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STUDIES OP EQUATORIAL AND POLAR LIGHT CHARGED PARTICLES 
EMITTED IN NUCLEAR EISSION AND DETERMINATION OP SCISSION 

CONPIGURATION 

A JIT KOMAR SINHA 
Ph.D. 

Department of Physics 
Indian Institute of Technology, Kanpur 
September 1981 

Nuclear fission is one of the important manifestations 
of large scale collective motion of nucleons involving a 
complicated interplay between the single particle and collective 
degrees of freedom. Inspite of nearly four decades of extensive 
investigations, an unified description of this process has not 
as yet emerged which is able to command general acceptance. 

The observed distributions in mass, charge, spin, kinetic 
energy and excitation energy etc. of the fission fragments 
depend sensitively upon the motion of the fissioning nucleus 
from the saddle point to the scission point. The initial 
conditions of such a motion are believed to be set at the 
saddle point. Various fission theories follow this motion under 
different assumptions regarding the coupling strengths among 
the particle and collective degrees of freedom. 

Studies in the light charged particle (LCP) accompanied 
fission (LCPP) bear a strong potential of probing the scission 
point configuration. In addition, it is believed that some 
characteristics of the LCPs such as-the dependence of their 
yield on the excitation energy of the compound nucleus may 
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reflect the saddle point effects. 

The LCPs observed in ICPF have been categorized as the 

equatorial LCPs and the polar 10 Ps, The former are emitted 

preferentially at right angles to the fission axis iwhile the 

latter are emitted . prefer enti ally along the fission axis. 

Equatorial LCPs are believed to be emitted from the neck 

region at the scission point while situation in the case of 

the polar LCPs is . ambiguous. Due to extremely low probability 

of polar emission, the experimental study, of this mode of 

fission is quite difficult. More experimental -information 

comparing the characteristics of the equatorial and the polar 

LCPs will be helpful for a clearer understanding of the 

mechanism of the. polar emission. 

The first part of the present thesis aims at studying 

some of the characteristics of the equatorial and polar LCPs 

235 

in the neutron induced fission of U upon a change in the 
excitation energy of the compound nucleus. The experiment was 

done for thermal and 600 keV neutron energies. Such a range 
of neutron energy was chosen in view of the earlier observa- 
tion of the variation in the LCP yield with neutron energy in 

235 

the keV neutron induced fission of U. The experiment was 
done by using a novel geometry providing a simultaneous detec- 
tion of the polar and equatorial POPs and permitting use of 
large area fission source without sacrificing the angular 
resolution. Such a simultaneous measurement of the equatorial 
and polar emissions provides a mutual comparison of the two 
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processes and it is expected to provide a clue to the mechanism 
of the polar emission. To the best of our knowledge, such a 
study for the polar LOPs has been done for the first time. 

The second part of the thesis concerns with the deter- 
mination of the scission point configuration by using the 
method of trajectory calculations. The calculations have been 
performed under point charge approximation utilizing some 
recently studied correlations in the alpha-particle accompanied 
fission. Improvements have been made in the method by incor- 
porating the nuclear forces and fragment deformations and 
effect of such improvements in the determination of scission 
point configuration has been investigated. Attempt has been 
made to see if we can explain the polar emission by incorpora- 
tion of the nuclear interaction. 

First chapter gives a short introduction to some 
important experimental and theoretical aspects of the fission 
process. Importance of an accurate knowledge of the scission 
point configuration has been brought out. 

Chapter II summarizes the work done on the 
equatorial and the polar light charged particles (LGP) emitted 
in nuclear fission, with emphasis on the polar LOPs, Various 

models proposed for the polar and the equatorial LOPs are 
described. Need of further experiments is pointed out and 
motivation of the present experiment is presented. 

The chapter III presents the details of the geometry 
used for the measurements of the polar and equatorial LGP 
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yields, the associated electronics and data recording arrange- 
ment , the data analysis, results and discussion. Monte 
Carlo method , used for calculating the detection efficiencies 
of the polar and equatorial IGPs and their possible intermixing 
in the process of identification has been described. The 
results show a variation in the yields of the polar and equa- 
torial LCPs when the neutron energy changes from the thermal 
value to the 600 keV. These results have been discussed in 
the perspective of emission mechanism of polar ICPs, 

Chapter IV introduces the method of trajectory 
calculations. The problems associated with the reconstructing 
of the scission point configuration by trajectory calculations 
are discussed and a survey of the earlier works has been 
presented. 

Chapter V starts with the trajectory calculations 
performed by using the point charge model. The method used 
in this work allows a free variation of all the scission 
point parameters except for some restrictions on the inter- 
fragment distance. The calculations were done for different 
ranges of values of the interfragment distance attempting to 
test the validity of different models regarding the saddle 
to scission point motion. The results of such calculations 
are presented and it has been shown that by a proper choice 
of the range of interfragment distance one can explain nearly 
all the experimental correlations in alpha-particle accompanied 
fission including some recent results on the angular distri- 
bution of the alpha-particle as a function' of the total 
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fragment kinetic energy. A relatively compact scission point 
configuration is favoured by these calculations. 

Incorporation of the nucloar interaction among the 
alpha-particle and fragments and an approximate inclusion of 
the deformation of the fragments have been done in the trajec- 
tory calculations to see the effects of the approximations 
made in the point charge model calculations, especially in 
view of the compact scission point configuration predicted by 
this model. All the initial -parameters describing the scission 
point configuration are permitted a free variation and selec- 
tion of the initial phase space is done by using the experi- 
mentally studied correlations in the asymptotic values of the 
dynamical variables. Attempt is made to provide a proper 
time scale in the calculations so that the origin coincides 
with the actual scission point. The results obtained are 
presented and they indicate significant changes in the scission 
point configuration. Due to the presence of nuclear inter- 
action, some of the alpha-particle trajectories are bent into 
the polar region. The results obtained for the polar alpha- 
particles are encouraging and show reasonable agreement with 
the experimental results in view of the sensitive dependence 
of the calculated characteristics of the polar alpha-particles 
to the parameters of the nuclear potential etc. The need for 
further improvements in the trajectory calculations, e.g. more 
realistic fragment deformations, their time dependence, the 
frictional forces in the nuclear interaction etc. are discussed. 



Lastly, the chapter VI contains a brief suimaary of 
the results of the experimental work and their interpretatio 
and the results of the trajectory calculations. 



CHAPTER I 


MCLEAR FISSION 


1 • 1 Experimental Aspects 

The nuclear fission is a type of nuclear reaction that 
results in the splitting of a nucleus into atleast two comparahl 
parts. It was first identified by Hahn and Strassnian^ after a 
series of radiochemical studies of neutron induced reaction in 
uranium. Nuclear fission is a very complex nuclear reaction. 
Inspite of very extensive theoretical and experimental inves- 
tigations, an unified picture of the fission process has not 
emerged. 

An isolated nucleus with a given excitation energy may 
decay via vfarious channels like y^^^ission, nuc3.eon or nucleon 
cluster emission, p-emission or fission. The fission mode of 
decay assume significant probability relative to the other 
modes for heavy nuclei or for highly excited nuclei, light 
nuclei, viz. gold, bismuth etc., are observed to decay signifi- 
cantly by the fission mode at high excitation energy while 
fission even with zero excitation energy (termed as spontaneous 
fission) occurs for heavy nuclei (A ^250). Spontar).eous 
fission enjoys an unique position towards understanding the 
fission phenomenon as in this, the initial state of the 
fissioning nucleus is a pure quantum state and known accurately 
A fissioning excited nucleus can be produced by several 
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methods, such as, by absorption of nucleons(p or n) , nuclei 

(^H, ^He etc.)^, photons^, by partial absorption like in a 

fi 7 

direct reaction (d,pf)° or by inelastic collision (a,a*f)'. 

Excitation may also be impai'ted by a fast chang'ing electro- 
magnetic field (Coulomb excitation) like in electron induced 
fission®or muon induced fission^. These different modes of 
formation of the excited nucleus essentially produce 
their characteristic population densities of different excited 
states in the nucleus, 

1 . 1 A Different Modes of Fission 

A nucleus, undergoing fission, most predominantly 
divides into two nuclei of comparable masses. This mode of 
fission of a nucleus is called the binary fission. Possibility 
of a division into three (or even more) comparable parts 
exists. Such a mode of fission is very rare and has not been 
studied extensively^*^ and we shall not be discussing it any 
further. Another interesting mode of fission has been studied 
in which a light charged particle (like H, H, ■'^H, He, He 
etc. ) is produced in addition to the two heavier nuclei. In 
the following we shall discuss some of the experimental chara- 
cteristics of the binary fission and the light charged particle 
accompanied fission. 

( i ) Binary Piss ion 


A fissioning nucleus undergoes large defoi'mations 
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before reaching the scission point, a stage in the fission 
process when division of the nucleus is complete with negligi- 
ble nuclear interaction remaining between the resulting parts. 
The nuclei formed at the scission point accelerate due to the 
mutual Coulomb repulsion and are called the primary fragments. 

The primary fragments are highly excited and neutron rich 

—19 

(Excitation Energy ^ 30 MeV) and within a short time ( 10 

sec), some neutrons are evaporated from them, which are known 

as prompt neutrons^^. Some of the characteristics of the 

12 

prompt neutrons (e.g. , the angulai* distribution with respect 
to the fragment direction) indicate that a part of them might 
have originated from the fissioning nucleus itself while it 
was on its way towards the split. 

The nuclei fomed after iiie emission of prompt neutrons 
are called fission fragments. These are unstable -towards p” 
emission and from stable isotopes, known as fission products, 
after a series of p” transformations accompanied by the 
emission of the so called delayed y-rays and the delayed 
neutrons. Do-oxcitation of the primary fragments is also 
possible by prompt Y~® 2 iission, however, it becomes competitive 
with prompt neutron emission only when the excitation energy 
approaches a value near the neutron binding energy. On an 
average about 9 of the excitation energy is disposed of 

by gamma ray emission. In addition, y-rays take aWay a large 
part of the fragment angular momentum as indicated by the 
large anisotropies in the gamma ray angular distribution 
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13 

with respect to the fragiflents , 

Sharing of the nucleons among the fission fragments is 
a topic of immense theoretical interest and has been studied 
rather extensively. liucleon distributions (viz., proton 
distribution, neutron distribution or their combined distri- 
bution) in nuclear fission has been studied by several methods 

15 

like radiochemical analysis , simultaneous velocity or 

16 18 

energy measurements of the two fragments ~ , study of the 

18 

X-rays and the y-^ays from the fragments , mass spectroscopy 
of the fragments^^, the A ]ii_s method using gas filled ioniza- 
tion chamber'^^ etc. At low energies, the most probable mode 
of division for nuclei with A > 229 is the asyumaetric mode in 
which the fragments possess unequal numbers of nucleons. 
Consequently the mass distribution shows two humps correspon- 
ding to the light and the heavy groups of the fragments. Por 
lighter nuclei (with Z< 83), mass distribution is single 
peaked indicating a dominant symmetric mode, Por intermediate 
nuclei (like .Ra, Ac) mass distribution is triple humped. A 
similai' behaviour with increasing preference for symmetric 
division is indicated in the mass distribution of very heavy 
nuclei like francium^^. However, as the excitation energy 
is increased, the syuimetric mode becomes relatively more and 
more prominent in all studied cases. The yield distributions 
of the isobars and the isotones produced in the fission process 
have been studied by several workers^^. Isobars t-Jith ‘even 
atomic number are found to bo relatively more abundant as 
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compared to those with odd atomic numbers. Distribution of 
mass has been studied more extensively as compared to the 
distribution of protons (charge distribution) and the distri- 
bution of neutrons. Attempts have been made to s'budy the 
proton and neutron distributions and look for systematic 
trends which might be reflected more clearly in their separate 
distributions rather than in their combined distribution, viz. , 

the mass distribution. 

22 

kinetic energy of the fission fragments lia.s been of 
great interest for understanding the fission process. The 
observed kinetic energy is contributed mainly by the potential 
energy of Coulomb interaction between the fragments just after 
the scission, except for a small contribution due to their 
initial kinetic energy at the scission point. ■ knowledge of 
the latter contribution plays an important role in the 
understanding of the fission phenomenon, 

A lot of work is devoted in studying the correlations 
among the various characteristics of the fission fragments 
(mass and energy), the prompt neutrons and the prompt y-^ays 
produced in the binary fission. 

In addition, angular distribution of the fragments 
relative to the direction of the particle inducing the fission 
has been studied for many nuclei. Such a study essentially 
provides experimental informations regarding the population of 
the excited levels in the fissioning nucleus and dynamical 
aspects of the fission process. 



( ii ) Lifiht Charged Particle Accompanied Mission (LQPF) 


27 

In lOPF , the nucleus splits into two nuclei of 
comparable masses and an additional light charged particle 
(LCP), e.g. p, d, “^He, ^He etc., '^'He being the most 

frequently emitted LCP. 

No drastic difference has been noted in the mass and 
the kinetic energy distributions of the fragments in the 
LCPP and binary fission. The observed difference is accoun- 
table by considerations of the mass and the kinetic energy 
taken away by the IGP. This result points to the essential 
similarity of the two processes. 

Angular distributions of the LCP’s with respect to the 
fragments, which is found to be sharply peaked nearly at 
right angles to the fission axis, quite definitively indicate 
that they are produced in the neck region nearly contempora- 
neous to the scission moment. Energy distribution of the- 
LCP together with its inter correlations with the ai\'';;ular 
distribution and the fragment mass and energy have been 
exploited by the method of trajectory calculations for 
finding o^it the scission point configuration in LOPF, which 
as indicated earlier, should closely resemble the scission 
configuration in binary fission. 

Studies of the prompt neutrons from the fragments 
in LCPF have indicated that the fragment excitation energy 
is less as compared to that in binary fission. In addition, 
the shift in mass distribution in LCPP from that in binary 
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fission corresponds to the light fragment side in near 

symmetric region and heavy fragment side in the asymmetric 

region. The primary fragments corresponding to these mass 

regions are highly deformed at the scission point as 

indicated by the average number of prompt neutrons emitted 

from them. This point may be taken to relate the origin 

of LCP- emission with the fragment of larger deformation. 

The TCP's discussed so far are the ones emitted 

preferentially along a direction nearly at right angles to 

tho fission axis and are called the equatorial TCP's. While 

studying the angular distribution of the TCP's, Piasecki 
24 

et al. discovered that the yield of the TCP per unit solid 

angle shows a rising tendency in the fragment direction 

contrary to a monotonous decrease as anticipated from the 

angular distribution of the equatorial TCP's. These TCP's, 

which show a preference towards emission along tho fission 

axis, are called the polar TCPs. The emission probability 

of the polar TCPs is found to be more towards tho light 

fragments than the heavy fragments. The fragment mass and 

kinetic energy distributions in the polar TCP accompanied 

25 

fission exhibit characteristics similar to that in 

the binary fission or equatorial TCPP. Experimental bharac- 
teristics of the polar TCPs, unlike the case of the equatorial 
TCPs, do not provide any unambiguous picture of tho origin 
of the polar TCPs. 
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^ ^ Th-eoJ^ies of Nuclear Fission 

Promptly after the discovery of nuclear fission, 

? fi 

Pleitner and Prisch suggested that the liquid drop model 

(LCtyi) of the nucleus provides a natural explanation for 

27 

it. Shortly afterward Bohr and Is/heeler provided a more 
quantitative description of the fission process using the 
liquid drop model. According to this model, the fissility 
parameter x = (Z /a) decides the stability of the nucleus 
towards shape vibrations. It was found that for values of 
X less than a critical value, there exists a potential 
energy barrier towards fission, which is higher for lower 
values of x. Calculations of the potential energy of 
deformation based on LEM indicated a value larger than that 
obtained from the experimental data for the barrier height. 
In addition, the model favours symmetric mode of fission in 
contradiction to the experimental results on the fission of 
the actinides. The incorporation of the shell effects in 
the simple liquid drop model was recognised as a possible 
remedy. Several experimental findings (i.e,, no shift in 
the position of the heavy fragment peah in the mass distri- 
butions for different actinide nuclei etc. ) indicated the 

presence of shell effects. 

2 8 

Strutinsky provided a framework for incorporation 
of the shell effects in the calculation of the potential 
energy of deformation. His method consists of adding to 
the liquid drop energy a deformation dependent shell 
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correction. Nilsson model was used to calculate the single 
particle energy levels in the deformed nucleus. A smeared 
level density was obtained by averaging the calculated 
single particle level densities over a large energy interval 
( ;^ 8 MeV) so as to wash out the shell effects. The differe- 
nce of the total energy of the nucleus obtained by the use 
of single particle level density from that calculated by 
using the smeared level density was used for the correction 
to the LM potential energy, 

A typical potential energy curve as a function of 
the deformation for an actinide nuclei is shown in Pig. 1.1. 
Shell corrected potential energy curve shows a double humped 
structure to be contrasted from the curve based upon the LDM 
also shown in the figure. Such a curve is obtained from a 
map of the potential energy of deformation calculated for a 
set of deformation parameters. The curve shown in the Pig. 

1.1 reflects the variation in the potential energy along a 
contour in the multidimensional space of the defoimation 
parameters, chosen in such a way that for all points on the 
contour there is a minima in the potential energy along the 
directions at right angle to the contour. A maxima thus 

obtained by climbing a hill along a trenched path, is called 

SI 

a saddle point and- as suggested by A. Bohr in his channel 
theory of fission plays a very significant role in the fission 
process. Angular distribution of the fission fragments in 
particle induced fission and the existence of sharp fission 
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resonances in slow neutron induced fission of odd-M nuclei 
have been successfully explained by the consideration of the 
saddle point quantum states. In such fission reactions, a 
large fraction of the excitation energy of the nucletis goes 
into the deformation energy at the saddle, leaving the nucleus 
cold. The level structure at the saddle point should resemble 
that of heavy nuclei with ground state deformation and thus 
are widely spaced at low excitation. This accounts for the 
observed sharp fission resonances. Angular distribution data 
are explained by considering the saddle wavefunction and 
assuming the conservation of the angular momentum about the 
nuclear symmetry axis during the motion from the saddle point 
to the scission point. Gorrelation of the odd parity states 
at saddle point with the observed mass asymmetry has also been 
suggested based on certain observations regarding the peak 
to valley ratio of the mass distribution in the koY-neutron 
induced fission of -^^Pu etc. . Evidence of the saddle 

point effects on the light charged particle characteristics 
in LOPE induced by the keV-neutrons has also been suggested'^ . 

Motion of the nucleus from the saddle point configura- 
tion to the scission stage is considered to play a crucial 
role in deciding various characteristics of the fission. As 
the nucleus proceeds saddle point onwards, the energy latent 
in the deformation becomes available to the other degrees of 
freedom corresponding to the collective motion and the particle 
motion. The coupling of the fission degree of freedom (which 



corresponds to a collective motion in which the imminent 
fragments move apart) to the other collective and particle 
degrees of freedom is a point of debate among the various fission 
models. The adiabatic model assumes this coupling to be very- 
weak and accordingly the released deformation energy sho-uld 
appear solely as the kinetic energy of the translation of the 
imminent fragments. On the other hand, the statistical model 
takes a strong coupling among all the degrees of freedom and 
assumes that the nucleus remains in a state of quasi- equilibrium 
throughout its descent till the scission point. Thus the frag- 
ments are supposed to possess comparatively small kinetic 

energies at the scission point. Some fission models (like 
34 

Worenberg' s"^ ) take an intermediate stand and assume the 
coupling to be strong enough to have statistical equilibri-um 
among the various collective degrees of freedom only. Norenberg 
further assumes that coupling between the collective and the 
particle degrees of motion becomes strong the moment the 
energy of the collective motion exceeds the energy needed for 
breaking a nudeon— pair. This model has met several successes 
and is being explored further. 

Predictions regarding the various characteristics, 
such as the mass and the kinetic energy distributions of the 
fragments are made by either following the motion of the 
nucleus from the saddle to the scission point by solving the 
relevant equations of motion (the dynamical treatment^^) or 
under the statistical model assumption by concentrating upon 
the scission point configuration (the statistical treatment^^ ) . 
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1 lh.e Dynamical Treatment 

In the dynamical treatment one considers the motion of 

the nucleus from the saddle point to the scission point by a 

suitable parameterisation of the shape of the nucleus and 

solves the relevant equations of motion to predict the shape 

of the nucleus at the scission point. The dynamical model 

assumes a complete statistical equilibrium in the distorting 

nucleus till the saddle point as this stage in the fission 

process is supposed to be veryslow ( % 10 Sec. ) in comparison 

—21 

to the characteristic nuclear time ( '){, 10” Sec.). At saddle 
point, the energy is distributed among the various dynamical 
variables as given by the laws of statistical mechanics, A 
Monte Carlo technique is used to calcxilate the diotributions 
in the values of the dynamical variables at the scission point 
by solving the equations of motions starting from the distri- 
butions at the saddle point. From these distributions, the 
mass, kinetic energy, excitation energy distributions can be 
calculated . 

A knowledge of the various inertial parameters and 

nuclear viscosities for different collective degrees of 

motion (which determine the coupling between the collective 

and the particle degrees of freedom) and their dependence 

upon the nuclear deformation and temperature as well as an 

accurate potential energy map over the deformation space is 

essential for a successful dynamical calculation. A lot of 

57 

preSiSnt research is devoted to achieve this information. 
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In addition, a versatile shape parameterisation is also 
required^®. 

"55 

Initial dynamical calculations with some simplying 
assumptions have not met much success and thus highlighted 
the need of more accurate knowledge of the various quantities 
used in the calculations. 

■57 

Recent calculations of the deformation energy curve 
for some actinide nuclei indicate that near the deformation 
corresponding to the second saddle, the nucleus becomes 
unstable towards mass asymmetric deformations. Thus a possi- 
bility of two nearly parallel barriers beyond the second 
saddle point has been suggested one leading to a preferential 
symmetric mass division while the other leading to a preferen- 
tial asymmetric mass division. This concept of two modes of 
the fission decay has received considerable attention lately 

and some experimental evidence in its support has been 
27 

suggested 

1 The Statistical Treatment 

Under the perview of the statistical model, Pong^^ has 
developed an approach for predicting various fission charac- 
teristics, At any instant, during the descent of the nucleus 
from the saddle to the scission point, the assumption of the 
statistical equilibrium permits a determination of the confi- 
guration of the nucleus at that instant. Thus the scission 
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point configuration is determined by statistical methods 
without necessitating the solution of the equations of motion. 
For such a calculation, an accurate information oi various 
nuclear properties like the energy level densities at various 
excitation energies and deformations for relevant nuclides, 
the ground state masses etc, is needed. 

Due to comparative simplification arising due to the 
assumption of complete statistical equilibrium, the statis- 
tical approach provides a possibility of predicting a large 
number of fission characteristics, such as, the asymmetry in 
mass distribution for actinide fission, dependence of the total 
fragment kinetic energy and the prompt neutron yields from 
the fragments upon their masses, spins of the fragments, 
characteristics of the LOPF etc. Predictions thus made are 
in fair agreement with the experimental results and the 
remaining disagreements are, according to Fong, possibly due 
to inaccuracies in various nuclear data used in the calcula- 
tions, However, inspite of a number of predictions of 
various fission characteristics in near agreement ^^ith experi- 
mental results (which provide an empirical support to the 
statistical approach), a conclusive fundamental justification 
for the statistical assumption is still awaited. 

As indicated earlier, the motion from the saddle to 
the scission point may not be entirely adiabatic or entirely 
non-adiabatic. In fact, an entirely adiabatic motion will 
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tend to cause a fast tearing of the nucleus, resulting in a 
rapidly changing nuclear potential which may bring out particle 
excitations. This will imply a coupling with the particle 
degrees of freedom and thus invalidate a pure adiabatic model. 
Similarly, the strong coupling assumption of the statistical 
model will permit a slow saddle to scission point motion. But 
in such a slow motion, nucleons get enough time to rearrange 
to the lowest energy states which leads to a lowered value for 
inertial parameters. This may cause an acceleration in the 
descent of the nucleus making the statistical model assumption 
less valid. Thus motion from saddle to scission point may 
seem to alternate between slow and fast stages and exact 
picture is not clear. 

An accurate information on the geometric and the 
dynamic configuration of the scission point seems crucial in 
ascertaining the validity of assumptions underlying different 
fission theories. 

The phenomenon of light charged particle emission has 
been identified as an unique tool for determining scission 
configuration. Their characteristics are also expected to 
depend upon the saddle point configuration. A theoretical 
and experimental analysis of various characteristics of the 
ICPs bears a strong potential for testing various fission 
theories . 

Among the ICBs , two different catagories naiaely 
equatorial and polar IGPs have been found, characterised by 
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their angles of emission relative to the fragment direction. 
Motivation of our work was to study the dependence of some 
characteristics of the two types of LOPs on the saddle point 
configuration and to look for a clue so as to understand the 
mechanism of the polar lOP emission. 

In addition, attempt has been made to reconstruct the 
scission point configuration by the method of the trajectoiy 
calculations. Possibility of explaining the polar emission 
by trajectory calculation is also explored. 
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CHAPTER II 


L ICHT GHARCEP PARTICLE A CCOMPANIED ElSSIOI i 


2 • 1 Equatorial Li,s:ht ChaTi^ed Particle Accompanied Fission 

The phenomenon of emisoion of light charged particle 
(specifically, a-particle) in nuclear fission was first 
observed by Alvarez in 1943 (quoted from Farwell et al.^). 
However, published account of this phenomenon was given by 
Tsien . Light charged particle accompanied fission (conven- 
tionally called as ternary fission) occurs once in about 
300-500 normal fission events. The direction of emission 
of the particle was found to be preferentially at right 
angles to the direction of motion of the fission fragments 
(the adjective 'equatorial' is added to differentiate from 
comparatively small but significant number of 'polar' light 
charged particles which are found to be emitted preferentially 
along the fragment direction). This fact led to the belief 
that these particles are emitted close to the scission point 
in between the two fragments. Thus, this particular mode of 
fission was recognised as a potential probe for investigating 
the scission point configuration. 

Besides the a-particlef, which are most dominantly 
emitted, emission of other light charged particles (LGP) in 
fission such as the isotopes of hydrogen, other isotopes of 



helium and some heavier isotopes (viz., those of lithium, 
berrylium, boron, carbon, nitrogen etc.) have been i-eported. 

As the a-particles form nearly 90 /. of all the LOrs, a- 
paj“ticle accompanied fission has been studied none exten- 
sively compared to the- fission with the emission of other 
LCPs. 

Relative probabilities of emissions of some of these 

252 

particles in the spontaneous fission of Of and thermal 

255 

neutron induced fission of TJ along >rith the values of 
mean and full width at half maximum of their energy and 
angular distributions are presented in Table 2,1, The 
angle of emission of ICP is found to depend upon its energy 
and the mass ratio and kinetic energy of the fraguients. 

Various investigations with multiparameter data acquisition 
systems have been carried out to study the energy - angle 
correlations in ternary fission. The average kinetic 
energy of the a-particles as a function of its angle of 

4 

emission displays a minimum around the most proDo.ble angle . 
The angular distribution of the a-particle becomes wider, 
with mild variations in the most probable angle of aaission, 

5 

as the kinetic energy of the a-particles increases , The 
angular distribution of the a-particles is also seen to depend 
upon the fragment kinetic energy, being broader for lower 
values of the kinetic energy'^. Mean angle of emission of 
a-particles with respect to the light fragment is found to 
decrease with the fragment mass ratio^. Average kinetic 
energy of the fragments decreases as a fiunction of the 
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a-particle energy almost linearly . In addition, several 
other correlations among energy-angle observables have been 
studied in the a-particle accompanied fission. ^Studies 
for other LCPs also indicate similar behaviour^. These data 
provide important input to the method of trajectory calcula- 
tions employed for obtaining the configuration of tJie fissi- 
oning nucleus at the instant of LOP emission. A comparison 
of the characteristics of the alpha-particle accomjpanied 
fission and binary fission indicates that the scission 
configurations of the nucleus corresponding to these two 
modes of fission resemble each other closely. Thus the 
trajectory calculations provide informations relevant to 
the study of binary fission. 

Kinetic energy and mass distributions of the fragments 

in the ternary fission with LGPs-j protons, tritons and a- 

8—10 

particles, have been measured . The distributions in 
the kinetic energy of the fragments in the a-particle 
accompanied fission and binary fission ai'e found to be nearly 
identical if an account is talcen of the kinetic energy taken 
away by the a-particle. The observed fragment mass distri-" 
bution in a-particle accompanied fission shows featvires 
very similar to those observed in the binary fission, Vvhen 
normalised to the same height, the lighter sides of both 
light and heavy peaks nearly overlap while their heavier 
sides are displaced relative to each other by 4 mass units. 

A similar trend is seen for the triton accompanied fission. 
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However, no unique information is derivable from these mass 
distributions as regards the dependence of the emission proba- 
bility of LOP upon the fragment mass due to absence of any 
information of the relative nucleon contributions from the 
two fragments to the LCP, This point, as discussed by 
Halpern^^, prohibits an assessment of the one step or the 
two nature of the LCP accompanied fission on the basis of the 
fragment mass distributions. The mass distribution in the 
proton accompanied fission shows some distinctive feature as 
compared to the a-particle accompanied fission. The observed 
distribution shows a relative preference for the symmetric 
division and both the light and heavy peaks ai‘e shifted 
towards the symmetric valley. Another dj^tinctive feature 
for the proton accompanied fission concerns with the compara- 
tively very broad angular distribution of protons which is 

interpreted to imply a surface emission mechanism for these 
1 2 

particles . If we assume that the surface emission is 
realistic, the observed preference for the ny.mmetric division 
appears qualitatively in the right direction as the symmetric 
division corresponds to comparatively more stretched scission 
configurationswi th more surface area for emission(as indica- 
ted by a dip in the total kinetic energy in. symmetric 
f ragmentatio n ) , 

The net energy release in the fo3rm of internal 
ezcitation and kinetic energy of the particles is about 2,5 
MeV less in a-particle accompanied fission as compared to 
that in the binary fission. The internal excitation energy 



is obtained from the average energy and nuiaber of prompt 
neutrons and y-pti-otons emitted from the fissioning nucleus 
and is a measure of the total * deformation -- excitation' ' 
energy of the fragments at the scission point. A comparison 
of the internal excitation energy as a function of mass of 
fission fragments in the a-particle accompo,nied fission and 
binary fission indicates that the emission of a-particle 
occurs at the expense of energy of the fr.agments with high 
deformation - excitation energy. 

The f;ill of 5t2,5 MeV in the total energy release in 
the a-particle accompanied fission in comparison to binary 
fission is less than that expected on the basis of the 
average separation energy of the a-particle of % 5 MeY. 

This point has been suggested to indicate a preponderance 
of fragments with even number of protons or neutrons in 
a-particle accompanied fission which results in the release 
of additional energy in the form of the pairing energy. The 
observation that the energy for the a-particle emission 
dominantly comes from the fragment which has hi^er defor- 
mation - excitation energy forms the basis for the belief that 
the required energy for the emission of a~p;irticle is 

imp, anted to it at the expense of the deformation energy of 

11 

the fragments. Halpern has suggested that this transfer 
of energy could occur due to a fast changing potential that 
will be seen by the a-particle due to a rapid snap and 
contraction of the protuberances of the nascent fragments. 
Rough calculations based on this sudden snap model of 



Halpern have met paz'tial success in explaining the relative 
yields of various light charged particles, however, no 
quantitative calculation on these lines exist. The pre- 
scission evaporation model proposed by ilamaiina et al.^"^ 
predict satisfactorily some characteristics of the light 
charged particles. However, a consideration of the energe- 
tics and the observed rather mild dependence of the light 
charged particle yield upon the excitation energy of the 
compound nucleus mates the evaporation mechanism as unten- 
able for explaining the light charged particle accompanied 

fission. An earlier suggestion of alp ha- emission from the 

15 

heavy fragment shortly after scission by feather has been 

discarded on the basis of the prompt neutron yield data of 
1 6 

Nardi et al. , as discussed in the review article by 
17 10 

Feather . Car j an has advanced some arguunents of quali- 
tative nature in support of pre-scission emission of the 
a-particle based upon rather elaborate calctilations perfor- 
med for the a-particle cluster-formation probability over 

19 

the surface of the nucleus near the scission stage . 

However, further calculations providing a quantitative 

support to the theory have not been done, 

20 

Fong has applied statistical theory for calcula- 
ting the probability of a-particle emission. He arrives 

at a value for emission probability of a-particle in thermal 

235 

neutron induced fission of U in agreement with the experi- 
mental value. The calculated probability for the emission 
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depends sensitively ujion the fragment excitation energy at 
the scission point. The observed decrease in the yield of 
prompt neutrons in a-particle accompanied fission as compared 
to binary fission is assumed to be entirely clue to a loss in 
the excitation energy of the fragments at scission, the 
deformation energy being assumed to be nearly same in the 
a-particle accompanied fission and binary fission. Apart 
from this assumption, the calculation suffers from the 
uncertainties associated with various parameters it requires 
for evaluating the emission probability. In addition, 
calculations have not been done in the case of other light 
charged particle accompanied fission. Such calcTilations may 
provide a better test for the statistical, model. 

Efforts have been made to study the dependence of 
the characteristics of the light charged particle accompa- 
nied fission upon the excitation energy of the fissioning 
nucleus. In the range of excitation energy upto ^ 20 MeV, 
a mild decrease in the yield of alpha-particle is observed 

while for still higher energy, the yield shows an increasing 
21-24 

trend' ' , The results of such experiments are difficult 

to interpret as there is an admixture of bao emissions from 

several nuclei at different excitation energies due to the 

occurence of second o.nd higher chance fissions. Several 

workers have performed experiments with excitation energy 

below the threshold of the second chance fission for a 

clearer picture of the excitation energy dependence of the 

25 

ternary fission probability. Results of uadkarni et al. 



for the neutron induced fission of with neutron energy 

varied from 1 MeV to 4 MeV in steps of 1 LeY showed that the 
a-particle yield in ternary fission is nearly constant. In 

2 fi 

a subsequent work, li'adkanni et al. found ttiat the yield 

of LGPs with kinetic energy above 12 he! jramains constant 

while the yield of LCFs with energy below 12 rieV varied 

significantly for the neutron energy region of 0.75 MeV to 

2 7 

1.75 MeV. Fluss et al. studied the yields of a~particle 

235 

and tritium from neutron induced ternary fission of U 

with neutron energy from thermal to 700 KeV. The alpha- 

particle yield was found to remain constant in contrast to 

a strong variation for the yield of tritons. Krishnarajulu 
2 8 

et al. .studied the energy distribution and yield of the 

LOP at several neutron energies in the range f 3 om thermal 

235 

to 1 MeV for the neutron induced ternary fission of U 
239 

and Pu. The alpha-particle yield showed a mild struc- 
ture in the neutron energy region around 200 KeV. They 
suggested that an explanation for the observed structure 
be sought in terms of properties (e.g. , pai'ity) of saddle 
point states of the .Missioning nucleus. Tlioy also presented 
results for the tritium yield but the errors were large due 

to the difficulty in proper identification in the absence 

29 

of the A E-S particle telescope. Sharma et al. carried 

out extensive* study of the yield of the protons, tritons 

and a-particles in the KeY neutron induced ternary fission 
235 

of U by employing a AE-E semiconductor telescope for 
particle identification. Their results confirm the findings 



2 R 

of K.rishnarajulu et al. , In addition, they found large 
increase in the yields of tritons and protons with neutron 
energy. 

2*2 Polar light Oharged Particle Accompanied fission 

One of the most important characteristic of the 

equatorial light charged particles (lOP) in ternary fission 

concerns with their sharply peahed angular distribution with 

respect to the fission axis. This observation provides a 

crucial argument for assuming that the IGP-emission occurs 

from the neck region nearly contemporaneously with the 

scission point of the fissioning nucleus. The Coulomb 

repulsion due to the fragments is expected to prohibit the 

LOP to emerge in directions close to the fragment direction. 

50-3 2 

However, several experiments indicated that there may be 

some IGPs produced in fission which are e.>flitted in direc- 

33 

tions very close to the fragment direction, Piasecki et al. 

carried out systematic investigation in the thermal neutron 

235 

induced fission of U for these close angle emissions. 

Their experiment established the existence of the so called 
polar ICPs emitted close to the fragment direction. The 
result indicated even a possible increase in the a-parttcle 
intensity (yield per unit solid angle) for directions in 
close vicinity of the fission axis. The mechanism of 
emission of the polar ICPs appejired as yet another problem 
before a fission physicist. As compared to the case of 
equatorial ICPs, the experimental study of the polar LGPs 



is quite difficult due to the lower probability of their 
emission (e.g. , polau' a-particle yield is about 1 /. of 
yield of the equatorial a-particles ) . The uncertainties 
associated with the measurements on the polci' LCPs as well 
as availability of onl.y limited data on them make it diffi- 
cult to assess the validity of various hypotheses advanced 
to interpret the polar LCP emission in fission. 

2,2a Experimental Studies in Polar Emission 

Observation of protons, tritons and a-particles in 

the polar direction has been reported^^"^^ for thermal neutron 

25 '^ 235 

induced fission of " U gjid U , spontaneous fission of 

^^^Cf and in the fission of by 42 HeV protons. 

However, no definite statement can be made regarding the 
dependence of characteristics of the polar LCPs upon the 
fissioning nuclei. In addition, no experimental information"^ 
existsso far on the dependence of the characteristics of 
the polar LCPs upon the excitation energy of the fissioning 
nucleus. Such data coiicoi'ning the characteristics of polar 
emission on parai'ietors of fissioning nucleu i could be very 
useful in understanding the mechanism of polar emission. 

Andreev et al. studied the yield dependence of polar 

q-particles for the thermal neutron induced fission of 233u and 
235 '(j. They report that ratio of probability of normal 
ternary fission and polar LOP accompanied fission are nearly 
same for the two cases, while the absolute ternary fission 
probabilities are quite different. 
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( i ) Anig^ilax Distribution, Intensity and xilnerA^y Distribution 
of Polar LCPs 

Several workers have attempted to obtain infoimation 

regarding the angular distribution of different polar LCPs. 

Due to the rarity of emission of polar LCPs, precise shape 

of their angular distribution is so far, not known. 

33 

Piasecki et al. carried out a study of the angular 

distribution of a-pai' tides over the entire range in the 

235 

thermal neutron induced fission of U, In the measure- 
ments for angles close to the fission a:8:is, the angle between 
lines joining the centre of the fissile source to the centres 
of a-particle and fission detectors was 0°, however, a 
significant probability existed for the detection of a- 
particles emitted within 20° from the fragiaent direction 
due to finite sizes of the detectors and the source. Angular 
detection efficiency curve was found to be ooaked at about 
12° with PVv'HL''I 16° . the angular distribution was recons- 
tructed from the experimental yields at di'Fferent angles 
by least squares method employing Monte Carlo technique. 

The reconstructed angular distribution (yield per unit 
solid angle) indicated a rise for directions in close 
vicinity (about 20°) of the fission axis, Adamov et al.^® 
investigated the angular distribution of tho polar LCPs 
using a semiconductor telescope for particle identification. 
Observed angular distribution showed a rise for extreme 
angles which was found to be larger for protons as compared 
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to that for a-particles, 

Recently Oaitucoli et studied .-uigular distri- 

bution of the a-particles in directions cloce to the fission 
axis. Their result indicates a mild increo,ae in the yield 
for directions along the fission axis. 

A. determination of the intensity of polar LGPs calls 

for an appropriate definition of the polar LCPs. Piasecki 

41 o 

et al. suggested a border-line of 25 for the angle of 

emission from fission axis for defining the polar LCPs based 

upon the shape of the observed angular distiribution. In 

addition, derivation of intensity of polar LCPs from the 

experimental, data demands a knowledge of their angular 

distributions due to finite experimental jaiigular resolution. 

Using the angular distribution data obtained in Ref. 33, 

41 

Piasecki et al. estimated polar a-particle intensity as 
5x10”^ of the equatorial a-particle in (n^^,f). Relative 

intensities of the different polar LCPs also need be correc- 
ted for the effects due to differences in thoir angular 
distributions. It is found that the intensity of polar LGPs 
emitted towards the light fragment direction is different 
from intensity towards the heavy fragment direction, being 
lesser by a factor of 3-4 for the latter case. These data, 
though yet to be precisely known, provide useful informations 
for understanding the nature of polar amission. The availa- 
ble data on the intensity of various polar LCPs is summarised 
in Table 2,2. 



Experimental intensity ratios for polar li^t charged particles 
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Energy distributions of the polar IGPs have been 
determined by various workers^^*”^^. The mean energies of 
the polar IGPs are found to be higher as compared to the 
corresponding equatorial LGPs while their dispersions seem 
to be comparatively smaller. In addition, the mean energies 
seem to be nearly independent of the fragment mass. Table 
2.3 lists the first two moments of their energy distribu- 
tions. A knowledge of the dependence of the polar LCP 
kinetic energy on its angle of emission could be useful for 
testing hypotheses regarding polar emission, but such data 
are not available. 

The kinetic energy of the fragments in polar LCP 
accompanied fission has been measured^^*^^. Average kinetic 
energy of fragments moving in the direction of polar a- 
particle is found to be less by about 15 HeV as compared 
to the corresponding value for the binary fission, while 
the decrease in the kinetic energy of the fragments moving 
opposite to polar IGPs is found to be much less ( 2 MeV), 

Part of the decrease in kinetic energy of fragments in the 
direction of polar IGPs is due to the recoil effect, while 
the change in kinetic energy of the oppositely moving 
fragments is expected to be comparatively uneffected by 
such recoil effects and should provide a good probe for the 
scission configuration in the polar LCP accompanied fission. 

Of interest is nearly linear fall in the kinetic 
energies of the polarly emitting fragments (that is, the 



Values of mean and standard deviation (in llel) of the energy distributions of polar 
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fragments moving in the direction of the polar LOP ) with 
increasing mass of the polar LCP (Hg. 2.1). This observa- 
tion seems to be relo,ted to a linear increase in the recoil 
with the mass of the polar LCP, which is expected on kinematic 
considerations. However, it is important to note that the 
effect of change of charge from to He on fragment recoil 
is somehow getting annul ed as there is not significant 
departure from the linear behaviour for this region. . 

Correlations between the kinetic energies of the 
polar ICPs and the fragments are not yet precisely deter- 
mined. However, average LCP energy is found to be nearly 
independent of the total fragment kinetic energy. Average 
fragment energy as a function of the LCP- energy is not 
reported so far. Such an information will be interesting 
in view of the observed anticorrelation in the equatorial 
LCP accompanied fission which plays an important role in 
ascertaining the scission configxiration. 

Total kinetic energy released in polar LCP accompanied 
fission (sum of the kinetic energies of the LCP and fragments) 
is not much different from the kinetic energy released in 
binary fission. In the cas6 of polar a-particle accompanied 
fission of U by thermal neutrons, the net kinetic energy 
released exceeds the corresponding value for binary fission 
by about 3 MeV"^ while in the case of polar proton and triton 

Obtained from Ref. 42 after averaging for polar emissions 

along light and heavy fragments. 
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accompanied fission, it is lesser ly about 2 heY. Thei^e 
results are similar to those obtained for the equatorial LOP 
accompanied fission. 

(ii ) Prafimont ilass Distribution in the Polar IQP Accompanied 
Fission 

Fragment mass distribution in polar ICP accompanied 

fission has been detoi'mined for the thermal neutron induced 

fission of and and spontaneous fission of 

252 

Of . Results on the thermal neutron induced fission of 
"^U by Piasecki et al.^ (and later confirmed by the work 
of Andreev et al. ) indicated a similaritj/’ in the mass 
distributions in tho polar and equatorial LGP accompanied 
fission. Average masses of the light fragments moving along 
and opposite to the direction of polar a-particles is found 
to be lowered by about 2.5 amu and 0.7 amu respectively as 
compared to the corresponding value in binary fission. The 
corresponding shift in the light fragment mass is about 2 amu 
for the equatorial a-particle accompanied fission. However, 
if we consider all the fission events for v/hich the a-particle 
is emitted along the fission axis within the polar cone, 
the average shift comes out to be 2.1 amu (obtained by a 
weighted average for tho emissions along light and heavy 
fragments). This number is in good agreement with the shift 
observed for equatorial emission. For tho proton accompa- 
nied fission of by thermal neutrons, a preference for 

the symmetric fission is revealed in both the polar emission 



as well as in the equatorial emission. 

Results of Nowicki et al.^^ for the x^olar a-particle 
accompanied spontaneous fission of ^^^Of indicate a very- 
small decrease in the mean mass of light fragments as compa- 
red to the corresponding value for the binary fission. This 

result is in disagreement with the corresponding result for 

252 

the equatorial a-particle accompanied fission of Gf. 

However, their result for the polar proton accompanied 

252 

spontaneous fission of Gf indicates an increase of about 
,3.1 am.u for the mean mass of light fragments (averaged for 
proton emission along and opposite to the fragment direc- 
tion) in the polar emission which may be compared with 

value of about 3.2 amu obtained for equatorial proton 

252 

accompanied spontaneous fission of Gf, 

2 • 2 B Hypotheses on the Hature of Polar i^lmission ' 

Unlike the case of equatorial LGPs, the characteris- 
tics of polar LGPs do not provide a clear indication of 
their origin and sever;'-! models to explain their emission 
have been suggested. Several alternate assi-imptions can be 
made as to the place and time of their origin in the fission 
process. For example, these LGPs might have been produced 
in the neck region at about the scission moment and later 
due to appropriate turning provided by the nuclear attraction 
of the fragments, they may emerge in the polar direction. 

Or, they might by diffracted into the polar region, as the 
quantum mechanical effects may play a significant role 



considering^ tho de Broglie wavelength of the a-particle at 
the scission moment. Another possibility is that these LCPs 
are omitted from the poles of the fragments (ends of the 
imminent fragments away from the nock) at time close to the 
scission, the omission being made possible due to sudden 
retraction of the nuclear matter at the poles. Yet another 
mechanism for tho emission of polar lOPs may be due to 
evaporation of the LCPs from the fragments. The evaporation 
of LCPs might bo significant if the fragment excitation 
energy is large, as for lower excitation energies neutron 
evaporation will overshadow the LOP emission due to the 
absence of a Coulomb barrier for neutrons. The evaporation 
being a compciTativoly slower process will occur from fully 
accelerated fragments. 

A largo number of hypotheses thus seem possible 
which can qualitatively explain the emergence of LCPs in 
the polar direction and have been discussed in fair detail 
by PiasGcki ot al."^^. 

One important question that seems relevant is whether 
the pol.ar LCP accompanied fission is a three particle 
breakup process or a two step process with LCi emission 
occuring from one of tho fragments after the binary division. 
In tho later type of mechanism, contribution of nucleons for 
the LCP would bo from the fragments moving along the direction 
of polar LCP (the possibility that polar LCP js emitted from 
the oppositely moving fragment is rejected a,s this would 
necessitate a very large kinetic energy of LCP at the 



instant of emission). On the other hand, for the first type 
of mechanism, the nucleon contributions from the fragments 
to the LGP can not be precisely ascertained if the LCP is 
emitted from the neck region. 

The only hypothesis which has been worked out in 
detail is the evaporation hypothesis'^^ of the polar 
emission. This theory has met partial success in explaining 
some of the characteristics of polar IGPs, such as the 
relative intensities of different polar ICPs and energy 
spectra of some of the LCPs (it fails in case of the energy 
spectra of tritons as well as of polar a-par tides emitted 
along the light fragments). The evaporation model permits 
a rather unambiguous prediction of the fragment mass distri- 
bution due to the assumed sequential nature of polar emission. 
Thus important test for the validity of evo,po ration hypothesis 
lies in the comparison of the predicted and experimental mass 
dio iributione for the polar IGP accompanied fission. lig. 2,2 
gives the required comparison. The observed disagreement 
seems difficult to be explained by the uncertainties associ- 
ated with the evaporation calculations. 

Apart from evaporation hypothesis, some calculations 

have been reported to test the validity of the hypotheses 

of nuclear orbitting of the LCPs in polar direction and 

quantum mechanical diffraction of LGP-wavepackets by the 

45 

fragments. Galculations by Dakowski et al. confirm that 
bending of the LGPs starting from the neck region due to 
nuclear attraction is possible with realistic values for 



Intensity (arbitrary unit 



ig, 2.2. Comparison of the measured (shown by open circle) ahd 
calculated (pre- emission) mass distributions of the frag- : 
merits moving along the polar LCP^ The bipartition datet ■’ 
are represented by dot-dashed line. The calculated curves 
are obtained by using various level density formulae emp- 
loyed in ; evc^orotjon model catculatiors : Truran et al C47 1 
(solid line) an^ Ignatyuk C403 (dashed tine). After 
-ail. ' ' 
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the dymimical coordinates of the fragments and the LCP. They 
calculated the nuclear force between the a-pexticle and 
fragments by using the a-nucleus optical po :;cntial (real 
part) and considered the effect of inclusion of phenomeno- 
logical friction term in calculating nucle'ir force. Fragment 
deformations at the scission point were considered. However, 
they conclude that a simultaneous reproduction of both the 
intensity and energy of the polar IGPs seems possible only 
if some dubious assumptions are made regarding the nuclear 
potential etc. CalculatioJis on the quantum mechanical 
diffraction of lOPs are reported by Andreev et al. and 

A fi 

hordyasz et al.^ . Kordyasz has done calculations by assuming 
the fragments as stationary and of equal mass. Calculations 
indicate that the observed intensity ratio ox the polar and 
the equatorial lOPs can be reproduced by such calculations. 

The calculated ratio of the polar and the equa.torial inten- 
sities is found to be strongly dependent upon the scission 
point configuration. As such these calculations appear 
quite promising, however, detrailed calculations with estima- 
tions of the effects due to various approximations made in 
solving the Schr'ddinger equation are not done. 

Lastly, we shall mention yet another model for the 

polar emission which assumes a delayed emission of the LOP 

( the delayed tripartion model described in the review paper 

of Piasecki et al."^^). Calculations have been performed to 

47 

test this hypothesis . These calculations with pure 
Coulomb force taken among the particles a,fter ICP-emission 
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have, however, failed to explain the polar emission. 

Thus we see that the present understanding of the 
phenomenon of polar emission is rather vague with a number 
of possible mechanisms prevailing. Experimental data on 
polar emission are limited. There is a need for more 
experimental data on the comparison of the characteristics 
of polar LCPs with those of comparatively well studied 
equatorial LCPs and on their dependence on the parameters 
of the fissioning nucleus. In the present work, we have 
attempted to study the dependence of the yields of polar 
LCPs emitted in the neutron induced fission on the neutron 
energy, Experiment has been done for the thermal and 0,6 
MeV neutron energies. The experiment is devised to provide 
a simultaneous measurement on the equatorial LCP yields. 

The possible comparison between the changes in the yields 
of equatorial and polar LCP with neutron energy assumes 
significance as the experiments on equatorial LCPs have 
indicated significant variations in their yields in this 
neutron energy region. 
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polar aud equatorial emission of LIO-HT oharged 

PARTICLES IN keV IISUTRON IRDIJCED FISSION 


3.1 The G-eometry 

The geometry used in the experiment is shown in Fig. 
3.1. The fission fragments were detected by a combination 
of two ionization chambers formed by division of a single 
ionization chamber into two parts by a multihole collimator. 
Such a geometry with two ionization chambers provides infor- 
mations regarding the direction of emission of the particles, 
which is necessary for a proper identification of the equa- 
torial or the polar nature of the detected LOP. The holes in 
the collimator form a hexagonal close packing structure and 
it is kept as close to the fission source as possible so as 
to minimise the uneffective area of the fission source. The 
separation of the nearest holes is chosen to ensure that the 
particles originating from a particular source area may pass 
through only one particular hole of the collimator. In 
addition, the collimator serves as a common electrode for the 
upper and the lower ionization chambers. The backing of the 
source fo.niis the collector electrode for the lower 
ionization chamber (LIC) while an aluminium foil serves this 
purpose for the upper ionization chamber (RIG). A ^ B-E 
semiconductor telescope for the detection of the light 
charged particles (LCP) is pl(aE^TBA)!viilS8AS^lminium 
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Fig. 3.1. Schematic diagram of the detection system. 
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foil. The aluminium foil is thick enough ( ^4 mg/cm^) to stop 
the fission fragments and the natural alpha-particles from the 
fission source. The chamber was filled with pure argon gas 
at a pressure of slightly more than the atmospheric pressure. 

Depending upon the point of origin and the direction 
of emission, a fission fragment trajectory may fall into one 
of the three general types indicated in Fig. 3.2. For clarity, 
only one collimator hole is displayed in the figure. For a 
trajectory of type (a) or (b), a pulse is obtained from the 
LIG while for (c) type trajectory pulses are obtained from 
both Lie as well as UIC. Possibility of a smaJ-1 charge 
collection by UIC for some trajectories of type (b) exists. 

The amount of this possible charge collection is minimised by 
increasing the distance of the aluminium foil from the colli- 
mator and applying a lower electric field gradient in the UIO 
as compared to that in the IIC. Most of the chad'ge produced 
inside the collimator hole due to trajectories of the type (b) 
either recombines due to very small electric field present or 
gets collected by the DIG. A consideration of the dimensions 
of the chamber implies a difference of atleast a factor of % 5 
between the DIG pulses originating due to the trajectories of 
type (c) and those due to the trajectories of type (b). In 
the experiment, the UIC piilses were found to be pretty large 
apart firom very small pulses due to above mentioned cause and 
the natural alpha-particles. These small pulses were cutoff 
by the electronic circuit used in the experiment. 



Al - foil 


Collimator 

Source 


(q) 



Fig. 3.2. Schematic illustration of typical fragment 
trajectories in the ionization: chamber. 





The pulses derived from the LIO, the UIG, the AB-S 
detector and the B-detector were looked for the coincidences 
by the electronic circuit described later. A three- fold 
coincidence between the A E-detector pulse, the E-detector 
pulse and the IIC pulse indicates that the fragment trajectory 
is of the type (a) or (b) and the LGP has passed through the 
collimator hole so as to be detected by the A 3-B detector 
telescope. Such a three-fold coincidence is taken to corres- 
pond to the detection of an equatorial LOP accompanied fission 
event. On the other hand a four- fold coincidence between the 
UIG pulse and the aforementioned three pulses indicates that 
both the LGP and the frag^ient have passed through the collima- 
tor hole. Thus such a four-fold coincidence is taken to 
correspond to the detection of a polar LGP accompanied fission 
event. 

3 . lA Choice of the Oollimator Size 

In deciding upon the size of a collimator-hole it is 
necessary to minimise possible contaminations to what we have 
designated as the equatorial events '(identified by the three- 
fold coincidence requirement) from the LGPs which are actually 
emitted in the polar direction. Similarly contribution to 
the LGPs which have been claimed by us to be due to polar 
emission (identified by the four-fold coincidence requirement) 
from the LGPs that are actually emitted* in the equatorial 
directions. Such contributions may be made arbitrarily small 




by chosing finer collimators. However, employing a very fine 
collimator means a reduction in the count rate resulting in 
near impossibility of doing such an experiment vrith the 
available neutron flux due to extremely rare emioaion of the 
polar LCPs. The dimensions of the collimator holes (Pig. 3.1) 
have been chosen, in the present experiment, to provide a 
reasonably low 'polar-equatorial' mixing keeping in view the 
need of maximising the 'vaJLue of the detection efficiency for 
the LCPs. A Monte Carlo calculation (described later) has 
been done to estimate the possible polar contamination of the 
equatorial events and vice versa. It is estimated that nearly 
7 , of the detected equatorial alpha-particles might be 

actually produced in the polar emissions. Corresponding 
number for the contamination of the equatorial protons by the 
polar protons comes out to be nearly 14 . Those numbers are 

the upper limit of the possible contaminations and have been 
calculated by using the equatorial and the polar data as 
given in the review paper by Piasecki et al.^. Similar 
calculations estimating the contaminations of the detected 
polar LCPs by the equatorial LCPs in our geometry indicate 
negligible contamination for the polar alpha-particles 
while the contamination of the polar protons by the equatorial 
protons comes out to be approximately 0.6 /. . Thus we see that 
with the present choice of the collimator size, the contamina- 
tions due to ' equate rial -polar mixing’ in our experiment is 
not significant (possible contamination due to scattering is 
estimated to be insignificant). 



57 


Figure 3.3 gives the relative detection efficiency for 
the LOP detected in coincidence with the UIO and the LIO pulses 
as a function of the angle between the LOP and the fragment. 

The efficiency is maximum around 14° and falls to negligible 
value after 34°. Thus the polar LGPs detected in the present 
experiment are emitted within a narrow cone a;iround the fission 
axis. 

The present geometry avoids the use of semiconductor 

fission detectors and permits the use of a broad fission 

source without sacrifying the angular resolution. Thus it is 

quite useful for studying the rare phenomenon of the polar 

emission in nuclear fission. In addition, a simultaneous 

study of the equatorial POPs is also possible wliich provides 

a direct comparison of the characteristics of the polar and 

the equatorial LOPs. A possibility of studying the angular 

distributions of the various LOPs exists by using collimators 

2 

of different hole-sizes. Work on this line is in progress . 

3 . 2 Experimental Details 

3 , 2A Details of the Double Ionization Chamber 

Voltage applied in the UIO and the LIO were adjusted - 
along with the time constants and the gains in the linear 
amplifiers to obtain nicely shaped pulses without reaching 
the saturation level of the amplifiers. The LIO pulses were 
found to be varying continuously from very low values to a 
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certain maximum value. Tlie low height pulses from the LIO 
contain contributions from the natural alpha-particles and 
the fragments emitted at small angles to the source and were 
not detected by adjusting the lower cutoff in the Timing- 
Single-Channel Analyzer (OETSG - 420). The UIO pulses were 
mostly large pulses apart from those due to the natural 
alpha-particles. In addition, the small pulses may also 
be, in part, due to the reason discussed in section 3.1. 

These pulses were discarded by suitably adjusting the lower 
cutoff value in the TSOA. The TSOAs were operated in the 
bipolar mode to obtain accurate timing signals from the LIC 
and the UIC pulses. With a known delay introduced in the 
UIC pulses, the fast negative timing signals from the two 
TSGAs were fed to the start and the stop inputs of a time 
to pulse height converter (TPHG, ORTEO-405) for studying the 
time relationship between the 110 and the UIC pulses. It was 
found that the two pulses (obtained when a fragment passes 
through the collimator hole) are coincident within 20 nS. 

The full width at one tenth of maximum of the peak in the 
spectrum of the TPHG output is nearly 80 nS. 

3 . 2B D etails of the LGP Detector 

A AB-B semiconductor telescope supplied by the 
ORTBG, Inc., USA, was used for the detection and identifi- 
cation of the LCPs. The AB-detector was a transmission 
mounted 49 p.m thick silicon detector with an active area of 


300 while the E-detector was 500 p.m thick wi'ch an active 
area of 450 The detector telescope was not kept very 

close to the collector foil of the UIO in view of the 
possibility of sparking due to the voltage applied to the 
collector. 

3. 20 Electronic Setup 

The electronic setup used for the signal processing 
and the data acquisition in the experiment is shown schema- 
tically in Pig. 3.4. The circuit was planned to store 
informations on the values of the AE and E once a LOP is 
detected coincident with a fission fragments. In addition, 
arrangement was made to provide information whether the 
detected fragment pulse came from the LIG alone or from both 
Lie and UIC. Pulses from the UIC, AS-detector and E-detector 
were fed into the respective low noise preampliiiers (ORTEO 
142A), which were kept close to the detection chamber. The 
preamplifier output pulses were transmitted to the linear 
amplifiers placed in the data acquisition room by shielded 
transmission cables of proper characteristic iiipedances to 
avoid reflection and distortion of the signals. The gains 
and the time constants of the linear amplifiers (ORTBG Models 
providing unipolar as well as bipolar outputs) were adjusted 
so as to provide nicely shaped pulses with their magnitudes 
within the working range (O V - 10 V) of the linear 
amplifiers. 
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The bipolar pulses of the linear ajnplifiers corres- 
ponding to the UIC and the IIG were fed to the timing single 
channel analyzers (ORTBO - 420) operated in the bipolar mode 
to achieve accurate timing signals. After a suitable delay 
in the IIC timing pulse to annul the time-mismatch arising 
due to a slight difference in the width of the IIO and the 
UIC pulses, the two timing signals (standard logic pulses 
of %0.5 pS width) are looked for coincidence by a slow 
coincidence unit (0E.TBC - 409). The output of the coincide- 
nce unit was algebraically added to the LIO timing signal by 
the ORTBO 'SUT-! and INVERT* unit (Model - 433A). The pulse so 
obtained will have a higher amplitude in case of a coincidence 
between the LIC and the UIC pulses and a lower amplitude 
otherwise. This two level logic pulse provided information 
whether the fragment passed through a collimator hole or not 
and Was recorded after inversion and proper sharping for each 
LCP-detection event. The AB-timing signal was -derived from 
the corresponding bipolar pulse by using the CROSSOVER 
PICKOBB UNIT (ORTEG - 407) which provides a fast-rising 
pulse at the instant of the zero crossing of the bipolar 
pulse. This timing signal provided the strobe input to the 
STROBED SINGLE CHANNEL ANALYSER (O.iTEC - 413). This arrange- 
ment provides accurate timing signal with adjustable lower 
cutoff value for the input pulses. The E-timing signal was 
obtained by using the CROSSOVER PICKOBB UNIT which 
itself provided the necessary lower cut off value ( 240 mV) 



for the input pulses. The pulse widths of the n and 
pulses were chosen such that the timing pulses derived from 
them were coincident with each other as well as with the 
logic pulse derived from the output of the 'SUM and INVERT' 
unit described above. These three pulses were looked for 
the coincidence by a slow coincidence unit (BOIL - Model 
LG666). The output of this coincidence unit generated the 
gate pulse to activate the data recording -unit. The unipolar 
pulses from the linear amplifiers corresponding to the AE 
and E detectors were brought in coincidence with this gate 
pulse by introducing the necessary time-delays with the help 
of the delay amplifier UNITS (ORTEC - 427 j EGIL - Da671). 

The data were recorded event by event on a magnetic 
tape with a three parameter data acquisition system which 
also provided a display of the 'channel numbers' of the 
three parameters. The three parameters were the A E-pulse, 
E-pulse and the two level logic pulse indicating the polar 
or the equatorial character of the detected LOP. 

3 . 2D T he Neutron Production 

Neutrons were produced using the proton beam from 
the Van de Graaff accelerator at Indian Institute of 
Technology, Kanpur, The accelerator is capable of providing 
a total of 100 pA current of protons with their energy 
adjustable over a continuous range extending up to a maximum 
of 2 MeV. The intrinsic energy stability of the machine was 




within + 5 heV and it could be improved further upto + 1 keY 
by using the slit stabilization s 3 '’stem fabricated in the 
accelerator lab. The machine was calibrated for the proton 
energy by using the sharp thresholds in the (p,n) nuclear 

7 

reactions with the li and the tritium targets. 

The Li(p,n)'Be reaction was used for generating 
thermal neutrons. The thermal! sation was achieved by using 
a 2" thick paraffin block. Thick lithium targets were 
made by vacuum deposition of the lithium metal over copper 
backings. After the evaporation, dry nitrogen gas was filled 
in the evaporation chamber so as to provide a protective 
lithiiom nitride layer on the lithium target. 

The ^H(p,n)^He reaction was used for generating 600 keY 
neutrons. The tritium target (obtained from the Isotope 
Division, BARG, Bombay) consisted of adsorbed tritium gas 
in the titanium layer deposited on copper discs, the ratio 
of the titanium atoms to the tritium atoms being 1:1. A 8 
curie tritium target was used which caused a spread of 
nearly + 160 keY in the energy of the neutrons. 

iTeutron target holder and the associa.ted beam 
transport system are shown in Figs. 3.5 and 3.6. Chilled 
water was circulated in the copper tubes covering the neutron 
target holder to avoid deterioration of the target due to 
heating by the proton beam. The aluminium grill arrangement 
(Fig. 3.6) at the exposure site reduced the possible conta- 
mination to the direct neutrons at the fission source from 
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Fig. 3.5. The neutron target holder, and collimator, . | 

I 




Aluminium Grill 



Fi§:3.6; Schemutic diagram of beam transport system. 
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the scattered neutrons. On the contrary, the chilled water 
cooling of the target holder may enhance such a contamination. 
However, overall thermal neutron contamination in the fast 
neutron exposure was monitored by the neutron activation 
technique using indium target and estimated to be less than 
0.01 •/. . 

The neutron flux was monitored by a standard long 
counter. The detector had a fairly flat response over the 
neutron energies above 1 MeY while its efficiency falls to 
nearly 70 */, of the maximum as the neutron energy reduces to 
the thermal value. However, the main purpose of the neutron 
counter was to monitor the deterioration, if a,ny, in the 
neutron target during the experiment. 

3 . 2jj) The Fissile Source 

The fissile source consisted of an uniform layer of 
enriched (93.7 /. ) deposited over an area of 4 cm^ on 

p 

a thick nickel backing. The source thickness was Jfe 5 mg/ cm 
and Was formed by using the spray technique in BARG, India. 

^ ^ E xperiment a l Procedure 

The iE and E detectors were tested for the energy 

241 227 

resolution by using thin alpha-sources of Am and Ac 
and the initial leakage currents were measured. The detec- 
tion chamber was tested for stability of the pressure of the 



filled argon gas thus ensuring the stahili'ty in the fission 
pulse height. The electronic circuit was tested by providing 
the pulses in various prepjnplifiers to simulate the polar and 
the equatorial LOP events. 

The A E and the E channels of the three 7oarameter 
data acquisition system were calibrated by using a precision 
pulse generator (ORTEG Model 204). The pulse generator was 
previously calibrated for the AE and E detectors. The pulse 
generator was also used to check the linearity of the electroni 
system over the complete range of the possible LOP energies and 
to get the lower energy cutoff values for the AS and E channels 
The stability of the electronic system was tested from time 
to time during the experiment with the help of the pulse 
generator. 

The lower level cutoff in the IIG and the UIC channels 
were adjusted to cut the pulses arising due to natural alpha- 
particles. The lower energy cutoff values were 0.7 MeV and 
1.5 MeV for the 'AE and the E channels respectively. The 
experiment was performed for the thermal and 600 keJ neutrons. 
The typical count rates ■( singles) of the IIC, the UIG, the 
A B and the E pulses are indicated in the Table 3.1. Typical 
equatorial LPG count rate was ^ 15/hour while the polar 
count rate was ^ 0.5/hr. These numbers refer to the thermal 
neutron run. Gount rate were less by a factor of 3-5 for 
the fast neutron run. Each fast neutron run was preceded 
and succeeded by a thermal run which provided a check on the 





possible drifts. 


1A3LE 3.1 

$ 

Typical count rates observed in the experiment for the thermal 

neutron run 


The Lie 

The UIC 

The E det. 

The E det. 

pulse 

pulse 

pulse 

pulse 

7600/sec. 

400/ sec. 

2.2/ sec. 

250/ sec. 


The leakage currents in the AE and B detectors were 
monitored regularly during the neutron exposure and the 
applied bias was adjusted to account for the fall in the 
voltage across the detectors resulting due to the increased 
leakage currents. This increase occurs because of the crystal 
defects introduced in the detectors due to the incident 
neutrons (damage being more severe for the fast neutrons) and 
causes a higher fall across the biasing resistance. The 
effect was more prominent for the B-detector due to its 
larger volume and so the biasing resistance in the preampli- 
fier (OETEG Model 142A) was reduced in this case. 

5 • 4 Data Analysis 

In the following, we shall be concerned with the 
details of the data analysis. The DEG-IO and IBM-ISOO 




computer of I.I.T, Kanpur were employed in the calculations. 


3 . 4A Particle Identification 

Particle identification (Pl) has been done by the PI 
signal generated from the A E and the B signals by using the 
following relation 

PI = (E + A E)^ - B^ 

with n being nearly 1.76. The E and ^E signals were derived 
from the corresponding channel numbers (stored event by event 
on a magnetic tape in the three parameter data acquisition 
system) and the calibration coefficients. 

The principle of the above method is based upon 
experimental range-energy relations for the H and He isotopes 
in the detector m e di urn ( that is, the silicon), which is 
summarised by the follO'Wing equation 

R(B) = B^ (1) 

where R(E) is range of the particle with energy E, n sfe 1.76, 

M = Atomic mass, Z = Atomic number and K is a constant (for 
energy E < 3 MeV, K is observed to deviate from the constant 
value). If the incident particle deposits an energy A B in 
the thin detector of thickness T and gets completely absorbed 
by the thick detector depositing an energy E in it, then we 
may write 
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T = (Range of the particle with energy E + A E) ~ (Range of 
the particle with energy E) 

or T K MZ^ [(S+AE)^~E^] 

Thus 

PI = (E+ aB)^ - E^ 

T 1 

= 7 . — 5 . 

MZ^ 

The above equation brings out the independence of the 
PI signal upon the energy of the incident particle. The 
importance of the method used lies in the fact that the 
need of a very thin AS detector is avoided which was an 
indespensable requirement in the earlier methods involving 
the measurement of B and (dE/dx). It is seen that the PI 
value depends upon the intrinsic characteristics (mass and 
charge) of the particle and thus is used for the particle 
identification, 

3,4B C hance Coincidence Contribution 

The present experiment is an extremely low counting 
rate experiment due to the rarity of the polar and the 

equatorial ICP emission. Thus care was taken to ensure a 
minimum of the contribution to the detected lOPs from those 
originating in nuclear reactions other than the fission 
reaction. This was done by imposing the three fold and four 
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fold coincidence conditions in the detection atri’angement . The cont- 
ribution owing to. chance coincidence has been estimated to be<l , 
of the observed number of the polar and equatorial lOPs. 

Hence no correction was made to account for the chance coinci- 
dence contribution. 

Detection Efficiency Calculation 

Calculation of the detection efficiencies of the 
equatorial and polar LOPs have been performed using the 
Monte Carlo method. Such a method employs the concept of the 
random number and simulates the physical events by chosing 
proper random variables. Various events, that' are to be 
simiuLated should be of statistical origin, i.e. , decay of 
the nuclei, absorption and scattering of the particles etc. 

In the present calculations, the physical events of 
the production of the LCP and the fragments were simulated 
mathematically in the computer for their motion aJfter the 
production. The direction of the motion were chosen by 
associating random variables having their probability 
distributions in accordance with the . corresponding experi- 
mental distributions. Among a large number of such 
'mathematical' events, in few the directions of the particles 
may be such that the conditions of the 'three fold' or the 
'four fold' coincidences are satisfied. The fraction of such 
events provides a good approximation to the detection effi- 
ciencies of the equatorial and the polar LGPs. In the 
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following, the procedure adopted is described in detail. 

1 • Simulation of the Light Charged Particle Accompanied 
Fission (IGPI) 

For simplicity, we shall consider a single collimator 
hole and the source area beneath it described by the shadow- 
cone (Fig. 3.7). 

( 3- ) Chosing the Point of Emission 

The point of occurence of the ICPF is chosen randomly 
within the shadow circle. For this purpose, we employ uniform 
unit randoni variable (to be denoted by RYU) which may have 
any value between 0 and 1 with uniform probability. Value of 
this random variable is supplied by the computei and we shall 
denote it by (RVU). If Z and Y are coordinates of the emission 
point, then 

X = Rj2iRYU) - 1] 

l' = R [2(RVD) - 1] 
s 

where R„ is radius of the shadow circle (Fig. 3.7) and (RVU) 

s 

used twice in the above equations denotes two independent 
values tah:en by the random variable RVU. 

Distance of the point of emission (= f(X^+T^)) from 
the centre’o'is calculated.. In case this distance turns 
out to be more than R„ , the current value is rejected and 

o 



prolxibrtity 
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new values for X and Y are chosen till the point of emission 
lies Within the shadow circle. 

The uniform choice of the point of emission is 
justified only if the fission source is uniformly thick and 
the neutron flux and energy do not change significantly over 
the source area. Both of these conditions are satisfactorily 
met in the present experiment. 

( b ) Oho Bing the Direction of Emission of the Fragments 

Direction of the fragment is specified by the spherical 
angle coordinates 0p. The coordinate system chosen (to 
be called later as space' fixed system) has its origin coinci- 
ding with the point of emission while the Z-axis is at right 
angles to the fission source. The experimental fragment 
angular distribution is isotropic in the thermal neutron 
induced fission while a small departure from isotropy occurs 
for the 600 keV neutron induced fission. However, this 
effect has been ignored and an isotropic emission of the 
fragments have been assumed. The angle Oj, has been chosen 
by using a random variable with its probability distribution 
proportional to account for the Sin © - dependence of the 
solid angle described by the directions (©,0), (©+d©,0), 
(©,0+d0) and (©+d©,d0). Tims random variable is given by 

4 

inverse cosine of an uniform unit random variable . The 
angle 0^, is chosen randomly within 0 and 2% with uniform 
probability. 
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( c ) Chosins the Direction of Emission of the L QP 

The direction of the LGP emission is correlated with 
the fragment direction. I'or accounting this, the direction 
of the ICP emission is specified by a proper choice of the 
spherical angle coordinates 0^ and 0^ defined in a coordinate 
frame with its Z-axis coinciding with the fragment direction 
(the origin coincides with the point of emission). 

The angular distribution of liie equatorial LOPs is 
peaJfced near 90° with PWHM (Pull Width at Half Maximum) 
varying from 20° to 90°. The polar ICP angulai’ distributions 
are not known accurately. They are found to be peaking along 
the fragment direction with PWHE4 of the order of 20° to 40°. 

The shape of the angular distribution has been taken 
to be represented by a Gaussian curve. Thus the angle 0^ is 
chosen with the help of a Gaussian random variable (to be 
denoted by RVG). The AVG is derived from the RVU by a 
procedure described in Ref. 4, figure 5.® shows a typical 
spectrum as generated by the RVG. In addition, a weight factor 
equal to Sin 0^ is also associated with each LOP event when 
the LCP-angular distribution is given in terms of yield per 
unit solid angle. The angle 0^ is chosen randomly within 
0 and 2tz with uniform probability. 

It is necessary to find the spherical angle 
coordinates of the LOP direction with respect to the space 
fixed system. This is achieved by using the rotational 
operator matrix. This matrix is obtained by taking the 



inverse of the transformation matrix for the Euler rotations^. 

After chosing the above coordinates, we calculate the 
coordinates of the following points : 

(i) Points where the fra,gment would have hit the lower and 
upper surfaces of the collimator. 

(ii) Points where the LGP would have hit the lower and the 
upper surfaces of the collimator and the plane of the 

AE and E-detector. 

Distances of these points from the axis of the 
collimator hole under consideration as well as from the axis 
of the A E and E detectors are calculated and the relevant 
distances are compared with the radii of the collimator hole 
and the A£ and E-detectors. Prom these comparisons it is 
discerned whether the LOP produced in the mathematical event 
is detected or not and if detected, would it have given rise 
to a three- fold (an equatorial event) or a four- fold (a polar 
event) coincidence. 

A large number of such mathematical events are 
generated by the computer, Number of the events which give 
rise to the detection of the equatorial or the polar lOPs 
are counted. If angular distribution used is given in yield 
per unit solid angle, a total ’weight’ is evaluated for 
various types (equatorial or polar etc.) of events by adding 
the ’Sin ©' weight factors. This number (or weight) when 
devided by the total n\imber (or weight) of the events gene- 
rated provides the relevant detection efficiency. 
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Calculations have been made on above lines by using 
the known angular distributions of the various LCPs. The 
experiment did not permit the identification of the light 
and the heavy fragment. In the efficiency calculations, 
this fact has been taken care of. The recoil of the 
fragments due to the LOP has also been included in the 
detection efficiency calculations. The detection efficiency 
depends upon the location of the collimator hole and the 
calculations have been done for different distances of the 
collimator hole from the common axis of the fission source 
and the LGP detectors. The distances of various holes in 
the collimator from this common axis are calculated and an 
average detection efficiency is evaluated. 

The thick fission source causes a loss of the energy 
for some of the fragments. This loss is severe for the 
fragments emitted in directions close to the surface of the 
source and some of these may not get detected in the ioniza- 
tion chamber due to the lower level cutoff set in the experiment 
Such an effect has been simulated in the calculations of the 
detection efficiencies by setting a lower angle limit (6^) 
for the emitted fragments. The value of the is chosen 
so as to provide a value of the detection efficiency for the 
equatorial alpha-par tides which gives the' equatorial 
alpha-particle yield in agreement with the accepted value. 



3.5 Results 


5 « 5A Particle Identification 

The particle identifier (PI) spectra for the LOPs in 
the thermal and the fast neutron induced fission are shown 
in the Pig. 3.9. A clear separation of the protons, tritons 
and alpha-particles is seen to occur in the PI spectra. 

The collimator used in the present experiment reduces 
the spread in the angle of incidence of the lOP on the AE-E 
telescope resulting in a better resolution in the particle 
identification spectrum. 

A simple Monte Carlo calculation of the PI spectra 
for the protons, tritons and alpha-particles has been 
performed for the geometry used. The angular spread of the 
LOPs, their energy losses in the gas and aluminium foil and 
the detector resolutions have been considered in the calcu- 
lation. The calculated widths of the PI peahs agree with the 
corresponding values in the observed PI spectra. 

The increase in proton yield for 600 keY neutron induced 
fission compared to the thermal fission is clearly seen in 
the PI spectra. 

3 , 5B Energy Distributions of the IGPs 

Figure 3.10 shows the energy distributions of the 
equatorial protons, tritons and alpha-particles correspon- 
ding to typical thermal and 600 keV runs. The energy 
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Fig. 3.10. Energy spectra of the equatorial light charged particles 
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distributions of the polar LOPs have not been included due to 
their poor statistics. 

The energy distributions for the ec^uatorial LOPs agree 
resonably with those obtained in earlier works (see table 2.1 
in Chapter II). However, a good quantitative comparison is 
not possible due to the low statistics in our experiment. 

The most probable energies of the equatorial protons and 
alpha-particles are nearly 8 MeV and 15 MeV respectively. 

A statistically significant value for the mean energy of the 
tritons could not be obtained. 

The proton energy spectrum indicates a sudden fall in 
the energy region near 10 MeV, This is because of insuffi- 
cient depletion depth of the E detector used in our experi- 
ment and thus for the protons incident on the telescope with 
energy more than 10 MeV, only partial absorption of the 
energy occurs. The alpha-particle energy distribution 
appears to be broader in comparison to the earlier results 
and indicates a somewhat double-humped structure. Such a 
shape of the energy spectrum is expected in our geometry 
due to (a) a lack of identification for the heavy and the 
light fragments and (b) discarding, of the fission fragments 
emitted at grazing angle to the fission source. The 
effect' arises essentially due to the angle-energy correla- 
tion for the alpha-particles. A Monte Carlo calculation 
of the alpha-particle energy distribution has been done 



82 


using the experimental results of Tsuji et al.^ on the values 
of the mean and standard deviation of the energy distribu- 
tion of the alpha-particles as a function of the angle of 
emission. Result of such a calculation is shown in Pig. 

3.11. The calculated energy distribution is found to have 
a structure with double humped character similar to that in 
the experimental spectrum. The widths of the calculated and 
experimental spectra also agree reasonably well. 

The energies of the polar lOPs were found to be 
consistent with the earlier results (see Table 2.3 in chapter 
II). The energy of most of the polar ICPs were found to be 
within twice the standard deviation from the accepted mean 
energies. 

3.50 Relative Yields of the LOPs 

Table 3.2 summarises the relative yields of the 
various LOPs. The yields of the equatorial and polar ICPs 
are normalised independently so as to give a value of 100 
for the corresponding thermal alpha-particle yields. The 
yields shown are after correction for the detection 
efficiencies of various IGPs. 

The detection efficiency is dependent upon the 
LGP angular distribution. Pig, 3.12 gives the calculated 
detection efficiency for the equatorial LOPs as a function 
of standard deviation of the angular distribution. 



Relative probability 






Table 5.2 


Relative yields of different equatorial and polar ICPs in the 

thermal and 600 keV runs 


Particle 

Equatorial 

Emission 

Polar 

Emission 

Thermal 

neutron 

fission 

600 keV 
neutron 
fission 

Thermal 

neutron 

fission 

600 keV 
neutron 
fission 

Proton 

1.2 + 0.4 

6.7 + .9 

5.4 + 5.5 

25 + 10 

Triton 

9.2 + 5.0 

10+5 

- 

- 

Alpha- 
parti cle 

100 + 9 

88+9 

100 + 53 

80 + 40 


The angular distribution data have been taken from the work 
of Tsuji et al. for alpha-particles while for protons and 

7 

tritons work of Raisbeck and Thomas was referred. The 
corresponding data for the fast neutron induced fission are 
not available. Hence, the detection efficiencies for the 
LCPs in the fast neutron induced fission have been calculated 
by taking the LCP angular distributions same as those in the 
thermal neutron induced fission. 

The angular distributions of the polar lOPs are not 
accurately known. In the detection efficiency calculations 
for them, the angular distributions were taken from the 
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references given in the review paper of Piasecki et al.^. The 
dependence of the detection efficiency of the polar lOPs on 
their angular distribution is shown in the Pig. 3.15. The 
efficiency shown is normalized with respect to the fragments 
detected in the upper ionization chamber (in the case of the 
equatorial LCPs, the efficiency was normalized with respect 
to the fragments detected in the lower ionization chamber). 

The statistical errors associated with the calculated effi- 
ciencies are of the order of 10 /. • 

The errors shown in the table 3. 2 are only statistical 
in nature. The possible errors arising due to the uncertain- 
ties in the detection efficiencies caused by the imprecise 
knowledge of the angular distributions are not included. 

The data on the polar tritons have been omitted due 
to their poor statistics. 

3 . 5D Comparison of Yields in Past and Thermal ITeutron Induced . 

Pission 

Table 3.3 contains the yields of the various IGPs in 
the fast neutron induced fission relative to the correspon- 
ding values for the thermal neutron induced fission. The 
errors indicated are only statistical as the errors due 
to uncertainties in the detection efficiency get cance- 
lled out in the relative yield calculations. 



87 


Table 3.5 

Yields of different equatorial and polar LCPs in the 600 keV 

run relative to that in the thermal run 


Particle 

Equatorial Emission 

Polar Emission 

Pro ton 

5.6 + 1.9 

4.6 + 3 

Triton 

1.1 +0.5 

- 

Alpha-particle 

0.9 + 0.1 

0 .8 +0.4 


Underlying such an observation is our assumption that the 
angular distributions do not change substantially when the 
neutron energy increases to 600 keY from the thermal value. 
Another possible factor which may affect the relative yields 
is duo to a change in th(3 angular distribution of the fission 
fragments with neutron energy. A calculation of the detection 
efficiencies indicates that due to this effect, a maximum of 
about 7 /. decrease might occur in the equatorial lOP detec- 
tion efficiencies while the polar LGP detection efficiencies 
remain nearly unchanged. 

It is observed that in fast fission, the yields of 
the equatorial and polar protons increases significantly 
while those of the equatorial and polar alpha-particles 
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remain nearly constant ¥ith. a slight decreasing tendency, 

Ihe yields of the equatorial tritons are nearly unchanged in 
the fast fission as compai’ed to the thermal fission case. 
Corresponding data for the polar emission are not presented 
due to their poor statistics. 

3 . 5-E The Absolute Yields 

ii'om the observed yields of the equatorial and polar 
alpha-particles in the thermal neutron induced fission, 
their absolute yields are derived by using the calculated 
detection efficiencies. These calculations are intended to 
provide a check for the experiment and a comparison with 
earlier works. The equatorial alpha-particle yield comes 
out to be 0.20 + 0.13 /. of the binary fission rate t , which 

O 

agrees reasonably with the experimental value of 0.17*/. •. 

Cross-section for the polar alpha-particle accompanied 

fission is evaluated to be 24 + 20 mb a number in reasonable 

2 9 

agreement with the result of the earlier work , Errors 
quoted include the contribution due to the uncertainties in 
in the detection efficiencies. The large errors in the 
evaluated cross-sections are mainly because of the uncertain- 
ties in the cut off value for the fission fragment pulses 

t 

This value is obtained chosing an appropriate value of the 
lower angle limit 'Oq' (see section 3.40) which also repro- 
duced the ratio of observed counts in the upper ionization 
chamber and the lower ionization chamber. 
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from the lower ionization chamber and the angulaf distributions. 
3.6 Discussion 

There is observed a manked increase in the yields of 
the polar and the equatorial protons for 600 keV neutron 
induced fission as compared to those for the thermal neutron 
induced fission (Table 3.3). The alpha-particle yield 
decreases slightly and the decrease is similar for the 
equatorial and polar emissions. The results on the equatorial 
alpha-particles and protons are similar to those obtained by 
Sharma et al.^*^* 

The most significant result of the present experiment 
is the observed similarity in the variations of the yield of 
protons and alpha-particles with neutron energy for the polar 
and the equatorial LOPs. The yield of the polar pimtons 
increases nearly five times for 600 keV neutron induced 
fission compared to the thermal neutron induced fission. An 
identical increase is observed for the equatorial protons. 
Similarly, the change in the yields for both the equatorial 
and polar alpha-particle is consistent. It may be emphasized 
that because of the simultaneous measurements of the polar 
and equatorial DOPs, these results are independent of noimali- 
zation uncertainties connected with the neutron flux and 
detection efficiencies etc. Moreover, contamination to the 
data from the particles with non-fission origin is very 
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little. Reactions with the fission neutrons may result in 
particles emitted in coincidence with the detection of the 
fission fragments and thus cannot be discarded by the coin- 
cidence requirements set in the experiment. However, 
contamination due to fission neutrons (or even fission 
Y-rays) is seen to be insignificant as (i) the relative 
yields of the different LCPs agree well with the earlier 
measurements and (ii) the average number of the fission 
neutrons or y-ra.ja per fission do not increase substantially 
nor do their energy distributions change significantly when 
neutron energy changes by a small amount of 0.6 MeV to 
account for the observed changes in the yields. 

Different models of the polar emission predict 
different dependences of the characteristics of the pol^ 
LOPs upon the excitation energy of the fissioning nucleus. 

In particular, nearly identical variations observed in the 
yields of the equatorial and the polar LOPs when the 
excitation ener'gy changes by a relatively small amount of 
0.6 MeV imply a similarity in the mechanisms for the produc- 
tion of the equatorial and the polar ICPs. The polar 
emission models, the nuclear orbitting model and the 
quantum mechanical diffraction model, possess this property 
and are thus favoured by the present experiment. 
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CHAPTER IV 


TRAJECTORY OALGUIjATIOTjS 


4 . 1 Introduction 

An experimental derivation of the coiifigu-ration^ 

3 . 

fissioning nucleus at the scission point (which sssentiald,^ 
means the kinetic energies of the fragments and their rel^^^ 
tive separation) is one of the important problems before 
fission physicist. If the fragments are found to have latge 
kinetic energies at the scission point, adiabatic model viQ_3_ 
be more appropriate in describing the fission process beyo^^^j 
the saddle point. On the other hand, a statis tical model 
will find support if the scission point icmetic energies 
out to be small. 

Results on the binary fission have failed to provi^jg 
a solution to this problem. Studies in the light charged 
prarticle accompanied fission (LOPE) indicate a Possibility 
of obtaining the scission point configuration by 
characteristics of the LOPE. This is possible as there are 
strong evidences to conclude that the LOP i® odncei at 
instant almost coinciding with the scissio^i point. Uso 
characteristics of the LGPF and binary fission indicate that 
their scission point configurations are very similar to e^ch 
other. The method based on classical mechanics to achieve 
end is known as the method of trajectory ^stions. 



The method of trajectory calculations essentially 
attempts to find a set of values of the dynamical coordinates 
( position— momentum ) characterising the scission point confi- 
guration (SPG) which when time-developed results in an 
asymptotic distribution of the dynamical variables in agree- 
ment with the experimental results. The method has provided 
very useful information regarding the SPG. However, various 
trajectory crilculations reported so far have not been able 
to provide an unambiguous picture of the SPG, Apart from 
the effects due to several simplifying assumptions made in 
c.alculation of the time development of the system, a forma- 
lism for reconstructing the SPG from the experimental results 
is not clearly laid yet. One of the most widely used method 
assumes a certain distribution in all or some of the SPG- 
variablos as specified by a certain number of adjustable 
parameters. These paxtameters are adjusted so that the 
calculated asymptotic distribution in the dynamical coordi- 
nates reproduces the experimental distributions. The distri- 
bution ch:-tracterised by the best fit values of the parameters 
is taken to correspond to the SPG. The method has obvious 
limitations due to use of a finite number of parameters in 
addition to the ambiguity associated with finding out the 
best set of values of the parameters as the search is made 
in a multi- dimensional space of the parameters. The second 
type of method attempts in some way to reconstruct the SPG 
directly by using the experimental distribution in the 
dynamical coordinates. The distribution in SPG-variables 



is not spocified a priori in these calculations. The method 
permits a search in the phase space of the dynamical coordi~ 
nates specifying the SPC without introducing any 'artificial' 
bias for any part of the phase space which may come up in 
calculations of the first type (mentioned above) due to a 
priori specification of the distribution in the SPG- variables, 

4.1 A Phase Space Gonsideration 

The experimental observables in the ICPP relevant to 
the trajectory calcifLations are the kinetic energies of the 
three particles and their relative directions of motion 
(which give their momentum coordinates). The SPG is essen- 
tially characterised by the initial momentum and position 
coordinates of the three p?ir tides, imong the experimental 
observables, it is the absence of knowledge of the position 
coordinates of the particles and the time elapsed from the 
scission point to the moment of the particle detection, 
which prohibits a calculation of the evolution of the three 
particle system backwrard in time to arrive at the SPG 
starting from the expei'imental distributions, 

Tho momenta of the particles attain constant 
asymptotic values in contrast to the position coordinates. 

The position coordinates, in fact, keep changing even 
asymptotically at a rate decided by the velocity of the 
particle. Thus, in the phase space (momentum-position space) 
'Of the three particles, points corresponding to various 
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configurations of the scission point keep translating with 
different characteristic speeds along the axes corresponding 
to the position coordinates. To present a raoie clear picture, 
we may replace the position coordinates (viz,, x coordinate 
of a particle) by a new set of coordinates obtained from the 
position coordinates by subtracting out a term consisting 
of a product of the time elapsed and the corresponding 
velocity components (viz., x-vt). The points in the phase 
space characterised by these modified position coordinates 
will now be stationary asymptotically. The coordinates 
'(x-vt)’ are independent of the position and time of the 
experimental detection of the particle and bear a signature 
of the dynamics at the beginning of the motion when the 
forces were significant. It is clear that an experimental 
measurement with window on the momenta of the particles deter- 
mines the integrated probability for various phase space 
points with same momentum coordinates but d.ifferent position 
coordinates '(x-vt)'. Thus the experimental data do not 
provide any means to distinguish between points in the phase 
space with same momenta but different '(x-vt)' coordinates 
and in the present calculation as a first approximation, we 
have associated equal probability weight to all such points. 
This probability is decided by the energies of the particles 
which in turn are related to their momenta. The expressions 
used for the probability weight are given in Chapter V. 

An additional problem of fundamental nature is associ- 
ated with the method of trajectory calculations. Given a 
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distribution of initial variables consistent with, the experi- 
mental data, wo may find another distribution of the initial 
variables consistent with the experimental data by a time 
evalutlon of the former. This leads to an ambiguity in the 
proper iduntification of the distribution of the initial 
variables with the scission point configuration or equiva- 
lently in finding the proper time scale with reference to 
the scission moment. 

4.2 Review of the Earlier Work 

1 

The pioneering work of Halpern brought out the 

relevance of the method of trajectory calculations for the 

determination of the scission point configuration (SPC). 

The possibility of using such an approach was, however, 

2 

already indicated by the work of Tsien , who attempted to 

investigate tho SPC using the experimental velocity distri- 

2'55 

bution of the a-particle in the ICPF of U. 

Till now, numerous attempts (atleast 24) have been 
made for reconstructing the SPC by the method of classical 
trajectory calculations. A large number of calculations 
start with a known distribution in some of the SPO-variables 
(while the others axe held fixed). As a result of the 
ambiguity associated with this type of approach and due to 
the general problem of chosing the proper time origin (as 
discussed above), several mutually contradicting results 
have been reported regarding t'he SPC. A part of the 



disagreement is also due to the errors in the earlier measu- 
rements on some of the characteristics (viz., the angular 
distribution of the a-particles) of the LOPP. i^ost of the 
calculations employ the experimental data on the a-panticle 
accompanied fission only, as this is the most common mode of 
the LGPP, In the following, a chronological description of 

various trajectory calculations is presented. 

2 

Calculations by Tsien were done with the initial 
momenta of the a-particle and the fragments taken as zero. 

He varied the fragment deformations and distance of the point 
of emission of the a-particle and could reproduce the experi- 
mental velocity distribution of the a-particle. 

Halpern^ worked out in detail the basic idea of 
trajectory calculations. He suggested that the use of 
classical mechanics in the time development of the three 
particle-system (a-particle and fragments) after the scission 
point could provide useful informations regarding the SPG 
and proposed the point charge model for the scissioning 
nucleus. His calculations indicated need of large kinetic 
energies of the a-particle (4.4 MeV) and the fragments at 
scission point for reproducing the experimental energy and 
angular distributions. He had taken a fixed value for the 
inter fragment distance (27 fm) and a fix:ed point of emission 
for the a— particle. Thus he considered a variation in the 
momentum coordinates of the particles but ignored that for 
the position coordinates in ascertaining the SPG. 
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b-eilikmon and iChlebnikov'^ performed trajectory calcula- 
tions to reproduce the energy and angular distribution of the 
a-particle under the framework of the droplet model. The 
magnitude of the initial momentum (whose direction was chosen 
isotropically) of the a-particle was varied for obtaining the 
best agreement of the calculated and experimental angular 
and energy distributions of the a-particle. The initial 
a-particle energy, uniformly distributed till a maximum value 
of 1 MeV, gave best agreement. However, the agreement was 
poor as far as the dispersion around the mean values (which 
were reproduced well) in the calculated distributions is 
concerned. 

Boneh er al.^, working under the point charge model, 

tried to reconstruct the scission point configuration by 

using the experimental data of Tfaenkel on the a-particle 
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accompanied spontaneous fission of Of. Based upon 
arguments concerning the angular distribution of the a- 
particle and the anticorrelation between the final energies 
of the a-particle and fragments, they prescribed high values 
for the initial energies of the a-particle (^ 3 heV) and the 
fragments Me?) at the instant of the emission of the 

a-particle. Under certain reasonable assumptions regarding 
the interdependence of various initial parameters characteri- 
zing the scission point configuration, they succeeded in 
reproducing the experimental angular and energy distributions 
of the a-particle. However, the arguments put forth by them 
were not conclusive. In addition, their conclusions would 
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have to be reconsidered in view of the disagreement between 
the result of Rraenkel^ and the presently accepted data on 

the a~pfirticle angular distribution, 

6 

ixatase calculated the a-particle trajectories under 
different initial conditions and obtained a parametric 
equation representing tho final a-particle energy as a fiinc- 
tion of initial parameters for different values of the mass 
ratio and total final kinetic energy of the fragments. He 
attempted to find the scission point configuration by 
employing the method of least square for reproducing the 
experimental a-particle energy distribution liJith fixed values 
of the final angle of emission of a-particle, mass ratio and 
energy of the fragments. The results indicated large values 
for the interfragment distance ( ^ 27 fm) and energies of the 
fragments 60 MeV). He used a certain functional form of 
distribution in the interfragment distance and initial 
a-particle energy. Apart from the approximations made in the 
calculations, the result may not be taken as conclusive 
because of the limitations arising due to a specific choice 
of the functional forms and possible ambiguities in the 
method of least squares. 

Strictly speaking, trajectory calculations aim at 
reconstructing the configuration of the scissioning nucleus 
at the instant of emission of the a-particle, which is 
expected to be very close to the actual scission point. This 
fact Was explicitly taken into account in the calculations 
of Raisbeck and Thomas'^. They assumed that the a-particle 
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emission occurs a time interval ’t^’ after the scission point 
(at which the fragments are formed with zero velocity). At 

the instant of eruission of the a-particle (with ’t ’ as 

—21 ^ 

0.4 X 10" sec), the average energies for the a-particle 

and the fragments were determined as % 2 MeV and 8 MeV 
respectively with inter fragment distance equal to 21.5 fm. 

The calculated energy distribution agreed with the experi- 
mental one while the calculated angular distribution was found 
to be narrower. In a later calculation by Haj o,gopalan and 

O 

Thomas , a simple modification in the initial position distri- 
bution of the a-particle was done to remove the discrepancy 
regarding the angular distribution. They extended their 

point charge model calculations for other light charged 

6 

particles, viz., pro tons, tritons and He and met partial 
success in explaining their characteristics wi'ch a single 
set of initial conditions, 

Blocki and Krogulski^ carried out computations for 
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the light chi-irged particle emission in the fission of H . 
They suggested that the experimental results are not suffi- 
cient to pin point whether the scission point configuration 
is compact or loose. Under the assumption of sudden approxi- 
mation, they considered the fragments to be at rest at the 
moment of scission. In addition to the a-particle, their 
calculations attempted to reproduce the energy distributions 
of the ^H, ^H, ^He and some heavier isotopes. In a 

subsequent publication^^ , they reported the results of their 
trajectory calculations for the light charged particles 
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emitted in bho spontaneous fission of these results 

included predictions regarding the energy and angular distri- 
butions of the particles. They concluded that the initial 
conditions for the emission of a-particle, tritons and ^He 
were quite similisr. However, an extremely elongated confi- 
guration was needed for the emission of protons and deuterons, 
which led them to conclude that the emission mechanism for 
these particles was different from that of a-particles, e.g. , 
emission from the surface of the fragments instead of from 
the neck. 

11 

Wardi et al. carried out trajectory calculations 
for the light charged particles ( n, ^He, “He and 

°He) emitted in spontaneous fission of ^ Gf, They attempted, 
for each of these light charged particles, to derive initial 
conditions that best fit the experimental results. The basic 
method of obtaining the initial conditions was similar to 
the one used by Boneh et al.,’^. The results indicated that 
both the initial distances between the fission fragments at 
scission and the initial kinetic energies of the particles 
tend to decrease as the mass of the light particle increases. 
The inter fragment distance in the proton accompanied fission 
Was found to be unusually large. They suggested that a 

different mechanism may he responsible for proton emission. 

12 

Trajectory calculations by Brtel could successfully 

reproduce the experimental results on the angular correlation 
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of the a— particle in the spontaneous fission of Gf. Values 
of % 0.5 MeV and % 1.0 MeV were obtained for the kinetic 



energy of the a-particle and fragments respectively. 

In the trajectory calculations ly Vitta^" and 
Fong , the initial conditions were taken directly from the 
predictions of the statistical theory of fission. They 

could reproduce some of the characteristics of the a-particles. 

15 

(iazit et al. used trajectory calculations to comple- 
ment their measurements on thermal neutron induced fission of 
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■'^U. The method adopted was similar to that of Boneh et al. . 

The calculations indicated a value of 24 fm for the inter- 
fragment distance with initial kinetic energies of the 
a-particle and the fragments as 3.4 MeV and 31 MeV 
respectively. , 

Musgrove"^ carried out trajectory calculations for 
different fixed values of the initial kinetic energy of the 
fragments and inter fragment distance. He permitted variations 
in the position and momentum coordinates in accordance with 
the uncertrdnty principle. Value of inter fr.agment distance 
equal to % 2 3.7 fm was found necessary for reproducing the 
experimental energy and angle correlations for the a-particle. 
Mean initial a-particlo kinetic energy was found to be ^ 2,75 
MeV. 

Tsuji et al.^'^ complemented their elaborate experiment 
on the energy and angular distribution of a-particles in the 
fission of by perfoiming trajectory calcxilations to 

derive the scission point configuration. A value of about 
22 fm for the interfragment distance was seen to be most 
consistent with their experimental results. The initial 



kinetic energies of the fr agents and a-particle were estima- 
ted to be ^ 10 ^’ieV and 1 MeV respectively. 

C:iLculations by Oavallari et al.^^, Possati and 
Pinelli ^ and Krishnarajulu and Mehta^*^ atterayt to achieve 
the initial conditions for the scission point configuration 
from the experimental measurements without any fitting or 
iterative procedure. 

18 

Oavallari et al. used a range of values for the 
initial p^iTameters, These ranges were divided into a number 
of intervals* The trajectories were calculated for a fixed 
mass ratio using different sets of initial values selected 
from these intervals and for five values of the initial 
fragment energy. Using the experimental a-particle energy 
distribution, they assigned a probability weight to each 
trajectory provided the final kinetic energies of the frag- 

4 

ments and final angle of a-particle emission fall within a 
specified most probable range. The initial conditions 
calculated by them indicate low values for the kinetic 
energies of the fragments and the a-particle. 

Possati and Pinelli^^ presented results of their 
extensive trajectory calculations based on the formulation 
presented earlier^®. The reconstructed scission point confi- 
guration for the most probable mode of fission showed a doupl 
peaked structure in the initial a-particle kinetic energy. 

The lower energy peak was situated near 0.5 MoV and corres- 
ponded to the prediction of the statistical theory of 


fission. 



Krishnarajulu and Hehta^® working under the assumption 
of point charge model, attempted to reconstruct the distribu- 
tions in the parameters characterising the scission point 
configuration by using the experimental kinetic energy distri- 
butions for the a-particle and the fragments. The input 
parameters took valj-ues distributed over a wide range. By 
using the statistical weight factor derived from the experi- 
mental energy distributions, they obtained the initial distri- 
butions, The Calculated, distributions, of initial energies of 
the fragments and the a-particle are broad and cover the 
predictions of both the statistical as well as the dynamical 
theories. 

21 

Gravron has done an extensive calculation to determine 
scission point parameters which fit experimental data on light 
particle emission from ^ Of, The parameters were given 
Gaussian distributions. In calculating the standard devia- 
tions, care was taken of the Heisenberg's uncertainty rela- 
tions. Trajectories were calculated for different sets of 
values of the initial parameters chosen randomly in accordance 
with their assumed distributions employing the Monte Carlo 
method. Apart from the energy and angle correlations for 
the a-particle, data of Oheifetz et al. on the emission of 
particle unstable ^He were reproduced satisfactorily, with 
interfragment distance fixed at 21,6 fm. The relevance of He 
data was brought out for getting over a basic ambiguity 
associated with the trajectory calculations which arises due 
to the fact that any set of scission point parameters can be 





extrapolated to later times to yield larger initial kinetic 

energies and distances. He extended his calcula,tions to 

explain the characteristics of other light particles, viz., 

1 3 

H and H. Triton accompanied fission could be fitted by 

using scission point parameters similar to those for the 

a-particle accompanied fission. Such a thing could not be 

achieved for the proton accompanied fission, 

23 

Ghoudhury et al. carried out the trajectory calcu- 
lations in a three point charge model for the a-particle 
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accompanied spontaneous fission of Of, By using the 
method of multivariate analysis and assuming the initial 
a-particle energy to be zero to find the earliest set of 
scission point parameters, they succeeded in reproducing 
various energy and angle correlations for the ternary fission 
of The scission configuration obtained by them 

corresponded to low kinetic energies of the fragments and 
a-particlo and a interfragment distance of ^ 20 fm. 

Trajectory calculations discussed so far do not 
consider the effect of the remnant nuclear interaction 
between the particles after scission. However, few of them 
have taken some account of this effect in so far as the 
a-particle was assumed to be absorbed by the fragments if 
the distance between them falls below a certain value. 
Incorporation of fragment deformation is also done in some 
of the works while calculating the Coulomb force. Recently, 
a critical survey of the point charge model of the scissioning 
nucleus is given by Dakowski et al. * They suggested that 
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the inclusion of the nuclear forces among the three particles 
in the ternjtry fission plays significant role in the deter- 
mination of scission point configuration. Their calculations 

21 

on the lines of G-avron indicate that due to inclusion of 
nuclear forces the interfragment distance is to he increased 
from 21,6 fm to 24 fm to reproduce the energy and angular 
distributions in the ternary fission. They emphasized a 
need of a more realistic model of the scissioning nucleus 
than the point charge model. 

25 

In a recent publication, Guet et al. have performed 

trajectory calculations with some account taken for the 

nuclear forces among the particles after the scission point 

using the experimental results in the thermal neutron induced 

26 

ternary fission of of Guet et al. . Their calculations 

indicated that the experimental correlations in the energy 

and angle of the a-particle and fragment energy can be 

successfully reproduced with a compact scission point 

configuration with interfragment distance around 21 fm and 

fragment kinetic energy ^8 MeV. 

Apart from the present review, reader is also 

27 28 

referred to the review articles by Feather , Asghar et al, 

29 

and ?andenbosch and Huizenga , 
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OHAPTER V 


THE PRESENT OUiOULATIONS 


5,1 Motivation. 

Importance of the light charged particle accompanied 
fission (LCPP) lies in the possibility of studying a large 
number of intercorrelations among the various characteristics 
of the LCPs and the fragments. These data are used by the 
method of trajectory calculations to reconstruct the scission 
point configuration (SPO) which plays a crucial xole in 
understanding of the fission process. The trajectory 
calculations attempt to find a set of values of the parame- 
ters specifying the SPO which can reproduce tne experimental 
intercorrelations between the various characteristics of the 
three particles. The configuration, specified by this set 
actually coriesponds to the instant of the light charged 
particle (LCP) production. The experimental characteristics 
indicate that the LCP is emitted at (or very close to) the 
scission point and thus the calculated configuration provides 
a good description of the scission configuration. 

The energy distributions of the LOP and the fragments 
and the angular distribution of the LCP relative to the 
fission axis have been studied for different values of the 
mass ratio for the fragments in the lOPF of several nuclei. 
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The LCP angular distribution is found to be peaked at nearly 
right angles to the fission axis (refer to Chapter II). It 
is, in fact, this important characteristic of the angular 
distribution which indicates that LCP emission coincides with 
the scission point. Several workers have studied the angular 
distribution of the LCPs, mostly concentrating- on the alpha- 
particles. Initial measurements^ indicated a large width of 
the angular distribution while the recent experiments with 
refined techniques prescribe an angular width (PWHM^l!18°) 
nearly half of the initial estimate. Trajectory calculations 
demonstrated the sensitive dependence of the calculated 
angular distribution upon the degree of compactness of the 
scission configuration. This led to more elaborate measure- 
ments. The angular distribution has been studied for some 
judiciously selected LOP events with the objective of narrowing, 
down the position-momentum spread of the SPO. Por example, 
the LOP-angular distribution has been studied by selecting 
ICPP events with the ICP energy restricted within narrow 
range^. Th(5 results indicated lesser focussing of the more 
energetic LCPs. Such a result has provided an additional 
clue to the SPO. Recently, a finer selection of the fission 
events has been made by Guet et al.^ by imposing the restri- 
ction on the fragment kinetic energy and studying the LCP 
angular distribution. They have also studied the LOP angular 
distribution with a window on the LOP-energy alongwith that 
on the fragment kinetic energy. They found that the mean 
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value oi the angular distribution of the alpha-particle remains 
nearly constant with an increase in the total fragment kinetic 
energy while the angular width shows a decrease with increa- 
sing fragment energy. The restriction on the fragment kinetih 
energy is expected to narrow down the variations in the initial 
position-momentum coordinates of the fragments. Keeping this 
in view, we have performed the conventional trajectory calcu- 
lations with particular attention to the angular distribution 
results of Guet et al. , 

In trajectory calculations, several approaches have 
been used in ascertaining the SPO (refer to Chapter IV also). 
Most of the present techniques start with a given distribu- 
tion in the values of the initial position-momentum coordinates 
of the three particles in LCPP. Such a distribution is charac- 
terised by a certain number of parameters and is assumed to 
correspond to the SPO. A distribution in the values of the 
final dynamical variables of the particles is obtained by 
solving the relevant equations of motion for different sets 
of values of the initial position-momentum coordinates. These 
sets are formed by a random selection of the initial coordi- 
nates so that the ensemble of such sets has a frequency 
distribution matching with the assumed distribution in the 
initial coordinates. The final distribution so calcxaated 
is brought in agreement with the experimental results by 
adjusting the values of the parameters characterising the 
initial distribution. The distribution characterised by 
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these adjusted parataeters is taken to correspond to the SPO 
of the fissioning, nucleus. The other mode of attack uses 
the experimentally determined distributions to provide a 
probability - weight for various combination of the values 
of initial dynamical variables possible for the SPC. The 
distributions in the initial values of the dynamical varia- 
bles (assumed to correspond to the scission point) are 
obtained with the help of these weights. • 

The later method appears more general and provides 
a search over the whole phase space possible for the position- 
momentum coordinates corresponding to the SPO. However, both 
of these methods suffer from the basic ambiguity regarding 
the fixation of the moment of the scission on the time scale 
chosen in the calculations. Assumption is normally made that 
the origin of the time scale corresponds to the scission 
moment. The wide and double humped distributions of the 

4 

initial kinetic energies obtained by krishnarajulu et al. 
and Possati et al.^ possibly reflect the effect of this 
ambiguity as pointed out by Pong^. In addition, several 
approximations are made in solving the equations of motion 
of the three particle, i.e. , the point charge approximation 
and ignoring nuclear forces among the particles after 
scission. Moreover, there is also a need of refining the 
overall method of trajectory calculations so as to arrive 
at an unique and possibly final description of the SPC, 





have made some efforts towards an improvement in the 
method of trajectory calculations. The calculations have been 

performed lor the alpha-particle accompanied spontaneous 
fission of 

5*2 P oint Qharge Model OalcuIa'tiQiis 

In this model, the fragments and the alpha-particles 
are substituted by point particles and their motion is followed 
classically under pure Goulomb interaction. A consideration 
of the de Broglie wavelengths of the particles involved indi- 
cates that they are small in comparison with the dimensions of 
the fissioning system. It implies that the classical 
dynamics should provide to a first approximation a good 
description of the motion of the three particles, 

5 • 2A Specifying the Initial Configuration 

for fixing the initial configuration of the three 
particles, a total of eighteen position-momentum coordinates 
need be specified. Here we ignore the internal coordinates 
of the particles as we are working within the framework of 
the point charge model. This implies a coverage of initial 
phase space region in an eighteen dimensional space. However, 
the specification of the positions of the three particles 
does not really require nine coordinates in view of the fact 
that we have to specify their relative positions only. ¥e 
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nsed only define three variables to achieve the specification 
of the relative positions, because the relative positions are 
fixed once the triangle formed by joining the particles is 
specified. Thus we need only three variables for an unique 
definition of the relative positions. In the present calcu- 
lations, the coordinates chosen are (i) two coordinates to 
specify the relative position of the alpha-particle and (ii) 
the distance between the fission fragments. Thus the number 
of the required coordinates gets reduced form 18 to 12. 

Purthermore, the assumption of planar motion implies 
that the components of momenta of the three particles at 
right angles to the plane through them are zero. Such an 
assumption may not be valid, for example, in case the ^alpha- 
particle has some rotational motion about the fission axis 
or the fragments have some rotational motion about the direc- 
tion of the alpha-particle emission, at the scission moment. 
The bending vibration of the fragments at right angles to 
the plane of the three particles may also invalidate this 
assumption. The first possibility can be discarded as this 
may imply a rotation about the symmetry axis (along which 
the imminent fragments separate) which is forbidden quantum- 
mechanically . The effects due to the rest of the possibili- 
ties are expected to be insignificant if the associated 
angular momenta are small. In addition, earlier calcula- 
tions^ have indicated that the assumption of the planar 
motion is a good approximation. Thus, we assume the planar 
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mo "b ion ’to be valid and the number of the required momentum- 
coordinates reduces to six from nine, The fissioning nucleus 
conserves its total linear momentum which is zero in the 
present case concerning with the spontaneous fission. Thus 
the total linear momenta of the three particles is zero which 
provides two equations of constraint for the momentum- 
components in the plane through them. This brings down the 
number of the required momentum variables to four. In the 
present calculations, we take two momentum variables to 
specify the velocity of the alpha-particle. The other two 
are to specify the velocity of the fragments. To achieve 
this, only the total kinetic energy of the two fragments is 
specified and the individual velocities are obtained by 
ensuring that the fragments have equal and opposite momenta. 
Such a procedure ignores the small perturbation due to the 
recoil of the fragments from the comparatively light alpha- 
particle. Thus the number of the required momentum- 
coordinates has been reduced to three. 

The parameters chosen to specify the scission point 
configuration are elaborated in Fig. 5.1 and described below; 

i) The coordinates X and T fix the position of the alpha- 
particle. The coordinate-system has its X-axis along 
the fission axis and the origin coincides with the point 
on the fission axis where the resultant electric field 
due to the two fragments vanishes. 
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ii ) Th.6 distrincB bstwoeix the fragpients# 

iii) The kinetic energy E° and the angle of emission 0° of 
the alpha-particle. 

iv) The total initial kinetic energy of the fragments (Ep. 

The Values of the various parameters mentioned above 

were varied over a wide range except for the coordinates Y 

and 6^. The coordinate Y was given a fixed value equal to 

zero while was chosen isotropically between 45° and 135°. 

4 7 

As earlier calculations * indicate that the final values of 
the dynamical variables are not affected significantly for 
a reasonable large change in the value of Y and because of 
limitations imposed on the computer time , such an assumption 
on the value of Y has been made. 

The range of the other coordinates were kept large to 
cover the values expected on the basis of the statistical 
and adiabatic assumptions regarding the saddle to scission 
motion. The values of the parameters are chosen so that any 
combination of their values has a total energy (kinetic + 
potential) within the (three sigma) range of the experimen- 
tally observed distribution in the total kinetic energy 
of the three particles. Table 5.1 gives the ranges of the 
initial parameters. 



Table 5.1 


The input-parameters 

and their ranges 

of values. 

Parameter 

Symbol 

Range of values 

Mass Ratio 

R 

1.0 to 1.8 

Interfragment distance 

B 

17 fm to 28 fm 

Initial distance of 
a-p article from zero 
field point 

I 

-7.0 fin to 6.5 fj 

Initial distance of 
a-particle from fission 
axis 

Y 

0.0 fm (fixed) 

Initial a-particle energy 

K 

0.4 to 4.0 MeY 

Initial fragment energy 


0.5 to 40 MeY 

Initial angle of emission 
of a-particle 

»1 

35° to 145° 


5.2b Mass and Charge Specification of the Fission Frafinients ' 

Calculations have been done for several values of the 
fragment mass ratio. The charges of the fragments are calcu- 

Q 

lated from their masses according to the equa.tion : 

Zjj = Mjj/ 2.587 and Zj^ = - Z^ 

where Z^,Zj^ represent the mean charges of the heavy and light 
fragments and Zrp corresponds to the charge of the fissioning 
nucleus minus the charge of the alpha— particle. Mg is the 
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mass of the heavy fragment. 

Etfect of the experimentally observed spread in the 
fragment charge is estimated to be small and has been ignored 
in the present calculations. The estimated sxjread in total 
potential energy due to this factor is less than 1 */. . The 
emission of prompt neutrons occurs after full acceleration 
of the fragments and its effect on the fragment kinetic 
energies is small. In the present calculations, the prompt 
neutron emission is ignored. 

5 . 20 The Forces Among the Particles and the Equations of 
Motion 

After specifying the initial position-momentum coordi- 
nates of the fragments and ICP, the classical equations of 
motion are solved to get the values of the final kinetic 
energies of the three particles and their relative directions. 

The forces, among the particles are taken to be purely 
electromagnetic, ignoring the effects of the nuclear forces 
acting in the initial stages of the post-scission motion. 

The dominant force among the particles is due to the 
Ooulomb interaction among them. The charge distribution in 
the fragment volume is not spherically symmetric as the 
fragments are deformed at the scission point. This implies 
that the electrostatic force contains contributions due to 
quadrupole and higher order interactions in addition to the 
dominant monopole interaction. Moreover, under the mutual 
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Coulomb tie.Ld, the fragments may get polarised and resulting 
dipole moment interaction may also influence the motion. 

However, their contribution to the total potential energy is 
estimated to be less than Z 4 Hence only monopwle 

interaction has been considered to keep the computation 
simple. 

The time taken for the full acceleration of the 
particles is of the order of Z 10 sec. : During this motion 
the electronic configuration in the atom re-uains unaffected 
and the fragments move in a stationary electron cloud. The 
electron cloud provides a force towards the centre of the 
atom on the fragments and the LOP. This effect is estimated 
to reduce the total energy of the three particlesonly by 0,02 */, 
and so it has been ignored, force due to magnetic interaction 
among the moving charged particles is smaller by a factor of 
(viV 2 /c^) as compared to the mutual Coulomb force. In the 
present c;aso, this factor is nearly 3x10 ^ and hence the 
magnetic force is not considered in the present calculations. 

The equations of motion to be solved are given by 


jyi ii 

^ dt^ 



i = 1, 3 

3 = 1 , 2 


where 



mass of the i"*^^ particle 

"t 

pQ 33 _-tion coordinate of the i particle 
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, ijli -f-v 

“• J coiaponent of the force acting on the i 
particle , 

These equations are integrated numerically on a 
computer with the help of following equations'^; 


At 


U- .• 

10 


(n+1) 


A t 


iT(n) 


U 


10 


(n) 

j,(n) 

M, 




At 


(n+l) 


X, 


( n+l ) 

ij 


U 


y(n) 


U. , 

10 


(n+l) 

iO 


p(n+l) 

v-^ 


A t 


At 


(n+l ) 


(n+l) 


where A t 


represents time interval 


a 



U 


ij 


is a constant 

represents component of the velocity of 

t h. 

the i particle 

t 

gives the mean velocity during the (n+l) 
step of the integration. 


The superscripts on ^ij correspond to the 

integration steps. 

The time interval has been chosen so as to be small 
in the initial stages when the forces are large and increase 
after each integration step. Such a choice saves computa- 


tion time 
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The forces are related to the position coordinates 
by the folio v/in/^ equation 


F. . 

13 


2 „ 


Ifcl 








t£{xa - 


The subscript k refers to the particles other than the 
one. 

The initial values (corresponding to the first step 
in the numerical integration) of the 1 . . and U. • are provided 

X J X J 

by the specified initial configuration. Initial time 
interval is chosen to be lO”^^ sec. while the 

constant ’a' is given a value equal to 0.1. 

The calculation was terminated after 81 time intervals, 
that is, after 3.3x10"’^^ sec., by which time the potential 
energy reduced to less than 1 /. of its initia-.! value. 

The final energies are given by 


h. 


1 

2 


M. 

1 


3=1 


and the an^^le the alpha-particle makes with the direction 
of light fragment given by inverse tangent of the 

ratio of the Y and X components of the velocity of the alpha- 
particle. 

A. simple consideration of the difference equations 
Used for integrating the equation of motion indicates that 
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alpiia-particle. The experimental distributions have been 
assumed to be of Gaussian shape. 

Distribution in a SPO-parameter is obtained by summing 
up the weights of all trajectories with specific values of 
the parameter. This procedure removes the possibility of 
any artificial, bias for any particular region of the initial 
phase space of the three particles. By using the statistical- 
weight factor alone a selection of the initial phase space 
is achieved which describes the SPC. 

( ii ) Restriction on the Interfragment Distance 

Different assumptions regarding the coupling coeffi~ 
cient among various degrees of freedom during the saddle-to- 
scission descent of the fissioning nucleus essentially 
predict different kinetic energies of the fragments at the 
scission point. This, in turn, implies different values of 
the inter-fragment distance because of the constraints 
introduced by the experimental kinetic energy distributions. 
Trajectory calculations have been performed to test the 
validity of various assumptions regarding the s addle- to- 
scission motion by restricting the interfragment distance 
D within a range of +0,5fto at values varying from 19.5 fm to 
26.5 fta. The range of variation of D chosen accounts for 
nearly 60 /. of the variation in the final kinetic energies 
of the fragments. All other SPC— parameters are allowed to 
have their normal statistical distributions as demanded by 
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the experimental kinetic energy distributions used for assig- 
ning weights to the trajectories, 

*;/e found that the trajectory calculations with all the 
SPC-parameters having their normal statistical distributions 
could not reproduce the angular distribution results of Guet 

3 

et al. . The experimental anticorrelation between the final 
kinetic energies of the fragments and alpha-particle could be 
reproduced^' with value of D restricted in the narrow range. 
These observations have further prompted us to make a detailed 
trajectory calculations with values of D restricted within 
a reasonable range around mean values lying between the 
limits set by the adiabatic and the statistical assumptions, 

5 , 2E Results and Discussions 

Figure 5.2 gives the angular distribution of the alpha- 
particle for mass ratio R equal to 1.4 with tvo windows on the 
total kinetic energy Bp of the fragments and permitting a free 

I 

variation in all the SPC-parameters. Mean value of the angle 
of emission of the alpha-particle with respect to the light 
fragment is unchanged for the two Bp windows while the 

B¥HM is more for the higher values of Bp. The result on 
agrees with the result of Guet et al. while the variation 
in the calculated BWHI4 with Bp is in opposite direction than 
observed in the experiment. In addition the calculated 
widths are higher than the experimental value of Bii/HM ^ 20 , 
Similar results have been obtained for other values of the 




Relative probability 
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Fig, 5.2. Angular distributions for Rs1,4 .with no restriction 
on th^ value of mterfragment distance (0) for 
different frogrnent eitergy,IEpl (Lines arc visual fits) 
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mass ratio. 

In our calculations, we have taken an isotropic distri- 
bution for the direction of emission of the alpha-particle. 

So the width of the angular distribution is expected to be 
more if Ep increases as then the Coulomb field responsible 
for the focussing of the alpha-particles diminishes at a 
faster rate. The effect of the interfragment distance on the 
angular width is essentially due to the subsequent limita- 
tions on the spread in the position X of the alpha-particle 
along the fission axis. 'At lower values of D, the spread 
in the position of the alpha-particle is expected to be less 
which should result in a reduced value for the angular width 
while the opposite effect might be expected for the higher 
values of D. The change in the calculated angular width 
with different value of Ep is mainly due to a combination of 
two opposite effects arising due to changes in the values 
of D and Bp. Higher Bp value corresponds to a reduced D and/ 
or an increased Bp. An increase in the calculated angular 
width with increasing value of Bp implies that the enhance- 
ment in the angular width due to increased Bp value is higher 
than the decrease expected on account of reduced D value. 

A possible reason for this may lie in the approximations 
made in the calculation. The selection of the spread in 
the X-value by the weight factor chosen may not be suffi- 
ciently realistic due to these approximations. It may be, 
for example, due to ignoring the nuclear interaction which 
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may cause more severe limitation on the spread, in the position 
of the alpha-particle (thus narrowing the angular* distribution 
more effectively) for the lower values of D. 

Krishnaradulu et al.^ found that other energy-angle 
inter correlations among the fragments and alpha-particle, 
viz,, variation of the with R, dependence of mean energy 
of the alpha-particle on etc, can be successfully 
reproduced without any restrictions on D, The experimental 
anticorrelation could be reproduced only after a restriction 
in the variation of D is imposed. 

We present results of the trajectory calculations 
perfoimed by fixing the interfragment distance at values lying 
between 19.5 fm and 26,5 fm in steps of 1.0 fm and allowing 
a variation of +0.5 fm around each value. Relations among 
the Bj,, for various values of R have been obtained 

for different values of D and compared with the experimental 
results. 

It is found after a detailed comparison of the 
calculated results with the experimental data that -only when 
the value of the interfragment distance is restricted below 
21 fm, the calculations reproduce all the experimental 
correlations. 

Tho calculations indicate that the observed narrowing 
of tho angular distribution with higher Bj, values and the 
anticorrelation between E^, and B^ provide a sensitive check 
for finding the correct value of the interfr agment distance. 



Relative probability 
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The angalah distribution of the alpha-particle for different 
values of the mass ratio has been found to be reproducible 
for more than one D-values studied. Similarly, the dependence 
of the upon could be reproduced for more than one 
D-value. Thus these correlations seem to be comparatively 
less effective in deciding upon an estimate of the actual 
inter fragment distance. 

Figure 5.3 shows the angular distribution of the alpha- 
particle for values of D restricted around 21.5 fm with two 
windows on the value o f Ep, The result for E restricted at 
21.5 fm are observed to be contrary to the experimental trend. 

3 

This indicates the significance of the results of Guet et al. 
in obtaining the SPO by the help of trajectory calculations. 

In the following we shall be describing various experi- 
mental correlations that have been reproduced by our trajectory 
calculations with D restricted below 21 fti. 

( i ) The Angular Distribution 

Figure 5.4 gives the angular distribution of the alpha- 
particles for R = 1.2 and R = 1.3. Similar results have been 
obtained for other values of the mass ratio. The calculated 
angular widths are in excellent agreement with the experimen- 
tal results^ 

Figure 5.5 shows the dependence of the mean angle 

of emission on the mass ratio. Also shovm are experi- 

1 10 

mental results of Fraenkel and Fluss et al. 


Our results 
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lie in between the two experimental results for as 3 mii 2 ietric 
mass division while for near symmetric region, the calculated 
Values are above the experimental result. This disagreement 
might be partly due to experimental inaccuracies associated 
with less probable symmetric fission. In addition, it might 
be a reflection of the approximation made regarding the 
deformation of the fragments, light fragment in near symme- 
tric fission is known to be highly deformed as seen by their 
enhanced prompt neutron yields. In our calculations without 
including the fragment deformations a value of ^92° is 
expected for near symmetric division due to the symmetry 
of the electrostatic field of the fragments (in absence of 
the fragment-recoil, it would be expected to be 90°). A 
deformed light fragment will exert lesser deviating force 
upon the alpha-particle as compared to the undeformed heavy 
fragment. This may cause a decrease in the value of from 
92° making the agreement with the experimental result better. 

( ii ) Angular Distribution with Window on the fragment Energy 

The calculated mean angle and the angular width 
for two windows on the fragment energy Ej, are presented in 
Table 5.2. Pig, 5.6 gives the angular distribution for the 
mass ratio 1.4 with two windows on the Bj,. The value of ©^-j^ 
remain nearly constant v/ith an increase in Bj, while the 
angular width shows a definite decrease with an increase in 
Ej,. These results are in good agreement with the experimental 






135 


results of Guet et al 


3 


Table 3.2 


Calculated angular distributions for various values of the 


mass ratio 

(R) and the kinetic energy of the 

fragments 

(El.)- 


165 + 

5 Me7 

175 + 

5 MeV 

R 

®al 

BWHM 

®aL 

BWHM 

1.2 

88.9° 

15° 

87.4° 

12° 

1.3 

85.5° 

17° 

85.3° 

14° 

1.4 

84.3° 

16° 

84.2° 

12° 


■5 

Guet et al. inferred that their experiment favours a 
SPC with low initial, kinetic energies. In such a case the 
experimental spread in the fragment kinetic energy Bj, will be 
primarily due to variations in the inter fragment distance D, 

In our calculations, the contributions to the spread in E^, 
are possible due to fluctuations in D as well as in Bp. Both 
of these variables are allowed a free variation over a reaso- 
nable range except for the restrictions arising due to the 
statistical weight factor. As discussed earlier, the change 
in the angular width with increasing Bp arises due to a combi- 
nation of a focussing effect caused by a possible reduction 
in D and a defocussing effect caused by a possible increase 
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in Ep. Our calculations could reproduce tiie result of G-uet 
et al. only with lower D values inspite of the fact that all 
values of D were checked in an identical fashion. This may 
he expected as the focussing effect associated ia.th D- 
vahiation becomes relatively dominant for lovrer D values as 
compared to the defocussing effect associated with Ep 
variation. Moreover, it must be pointed out that the nece- 
ssity of taking low I)-value for reproducing the results of 
Guet et al. may become less critical once we remove the 
approximations of ignoring the fragment deformations and the 
nuclear interaction. 

( iii ) Average Energy of the Alpha-particle (B^) as a Function 

of its Angle of Emission ) 

Figure 5.7 shows the variation of the average energy 
of the alpha-particle as a function of its angle of emission. 
Also plotted are the experimental results of Tsuji et al.^^. 
The position of the minima in the calculated and the expeii.- 
mental curves are in good agreement with each other. The 
calculated and experimental values agree within 0.5 MeY 
for different values of 

(iv) Dependence of on Bp and Bp on E^ 

Figure 5.8 shows the correlations among the fragment 
energy Bp and alpha-particle energy B^ for the mass ratio of 
1.3. A linear relation is observed between average fragment 
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Fig. 5.7. Average energy of the alpha-particle Fee as a 
function of its ongle of, emission 9<3ei.. Also 

shown are the experimental results of Tsuji etal 
cm (Lines are visual fits). 
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energy E^, and (Big, 5.8a). The slope of this line comes 

out to be -0,4 in excellent agreement with the experimental 

12 

value of -0.44 obtained by Mehta et al. . The relation 

between average alpha-paJ’ticle energy E^ and -turns out 

to be nonlinear (Fig. 5.8b) in agreement with the experime- 
12 

ntal result . 

These correlations have been calculated for different 
values of R. Figure 5.9 gives the calculated anti correla- 
tion ( 3 Bp/ 3B^) as a function of the mass ratio R, Also 

12 

shown are the experimental results of Mehta et al. . The 
agreement between the two is seen to be satisfactory. Boneh 

7 - 

et al. suggested that the anti correlation between Ip and 

B^ can be reproduced only by chosing high initial kinetic 

energies 20 /, of their final values). On the other 

13 

hand Ghaudhury et al. could reproduce the anti correlation 
condition with smaller values of the initial kinetic 

4 

energies. The calculations of Erishnaradnlu et al. covering 
lower as well as hi^er values for the initial kinetic 
energies indicated that the correct anti correlation is 
reproducible if the D-variation is limited in a narrow 

4 7 

range. Necessity of limiting the D-variation is attributed * 
to a possible positive correlation between E° and Ep arising 
with D-variation which gives rise to a positive correlation 
between B and E™. In the present calculations, the 
permitted variation in D is reasonably wide accounting for 
nearly 60 /. of the observed spread in Ep. With this much 
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Fig. 5.8. ( q ) Dependence of average fragment energy Ep 
. on the' alpha -particle energy ; 

(b) Dependence of the average alpha -par ticle . 
energy Eoc cbi the frogment : energy Ep . Mass 
ratio is 13. (lines are visyal\fits} / V ^ 
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of vaxiafion in the D-value, the anticorrelation condition 
could not be reproduced for higher D-values. At lower D- 
values, the interplay between the various SPG-parameters 
in deciding the asymptotic energies and angles becomes 
fairly strong. Ttiis mahes any qualitative reasoning for 
searching the source of the observed anti cor relation rather 
difficult and less reliable and emphasizes the need of an 
accurate and unbiased approach for ascertaining the SPG. 

(v) Distributions in the fragment Kinetic Energies (£p) and 
the Ain ha-par tide Energy (E ^) 

Iwe used experimental Bj, and distributions for 
assigning the weights to -the trajectories. As a cross 
check for the procedure used, we have calculated the final 
fragment energy and a-particle energy distributions. Mean 
values of Ej, and are quite accurately reprodu oed for all 
the Values of mass ratio. The width of the Ej, distribution 
agreed quite satisfactorily with the one used in assigning 
the weights. However, the width of the alpha-particle energy 
distribution was low (about 0.6 times the width used in 
assigning the weights). This might be a reflection of the 
approximations involved with the point charge model 
calculations . 

( vi ) Distributions in the SPG-parameters 

figure 5.10 (a) and (b) give the distributions in the 
initial a-particle energy and the total fragment kinetic 



5,10. Cg) and, (b) art; the calculated dtstd^tions in the inlffal 
cj^" particle energy IE 4^} and fragnnent Jcinetic fhergy CEf ) 
, . ic) and Cd) are the distributions in E^ and Ep as .^tdinc 
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energy obtained from our calculations. The calculated distri- 
butions of the and are sharply peaked at about 0.5 MeV 
and 1.2 MeV respectively, wLth high energy tails. They seem 
to resemble very closely the Maxwellian distributions shown 
in Pig. 5.10(c) and (d). They have been obtained by assuming 
a Maxwellian energy distributions for the a-particle, light 
fragment and the heavy fragment with temperature 0.5 MeV. 
Distribution in B^, the sim of the energies of the two 
fragments, is obtained from individual distributions by using 
the Monte Carlo method. Such a shape of the initial kinetic 
energy distribution for the light charged particle was used 
by Block! et al.^^ and might be expected under statistical 
equilibrium. 

The distribution in the position of the alpha-particle 
(the X-coordinate ) is shown in Pig. 5.11* It is found to be 
peaked at about 1 fm towards the light fragment from the zero 
field point with a fall width at half-maximum (PV/HM) of about 
5 fm. This distribution actually represents distribution of 
the projection of the actual emission point on the fission 
axis. It may be inferred that this spread in the X-coordinate 
contributed most dominantly to the angular width in the 
present calculations. 

5 . 2P Conclusio ns 

The experimental correlations between the energies 
of the three particles and the angle of alpha-particle 
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X in fno 

Fig. 5.11. Distrbution of the point of emission CX) 

of the alpha -particle. The point with Xo=0 fm 
corresponds to the point of zero coulomb 
fields of the fragments; Une is a visual fit.> 
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emission could be successfully reproduced by our trajectory 
calculations. Within the errors due to approximations made 
in the point charge model, a compact scission point configu- 
ration with low initial energies is favoured by our calcula- 
tions. However, as pointed out earlier, there is a basic 
ambiguity in determining the SPG using trajectory calcula- 
tions. During the time evolution of the trajectory, any 
time along the trajectory can be taken as the scission point. 
However, we find that for a typical trajectory characterised 
by = 0.5 MeV, = 1.2 MeV, = 0 fm and D = 20.5 fm 
(these are the most probable values as determined by our 
calculation), distance of the alpha-particle from the 
fission axis, exceeds 3 fia when the fragment separation 
becomes more than 23 fm. Davies et al. give an estimate 
of the neck radius of about 2 fm at scission, Takiing 3 fsi 
as a limiting value of the neck radius at scission, we see 
that for D >23 fm, the configuration of the three particles 
is nonphysical for the scission point (assuming that the 
a-particle emission coincides with the scission and occurs 
in the neck region). Considering the fact that the scission 
point parameters can be extrapolated to later times which 
will yield larger values for the initial kinetic energies 
and the fragment separation, we can get estimates of the 
range of variations in the initial parameters by imposing 
the condition that Y at scission point should be < 3 fia. We 
find that under this restriction, upper limits of the possible 
values of E°, and D would be 1.4 MeV, 16 heV and 23 fm 
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respectiYely, 

5*3 ImproYements in tlie Point Charge Model 

¥e haYe, so far, discussed tiie results of the traje- 
ctory calculations ly using the point charge model for 
describing the post scission motion of the fragments and 
alpha-particle. ¥e could reproduce various energy-angle 
correlations among the three particles when we restricted the 
interfragment distance (D) around 20.5 fh. 

A simple consideration of the fragment sizes and 
range of the nuclear force indicates that the effect due to 
the nuclear forces may he significant for the alpha-particle 
emission points lying close to the fragments. Ihe Z-distri- 
hution obtained by the point charge model calculations (iig* 
5*11) shows a spread of nearly 5 fti. Ihis spread plays a 
significant role in reproducing the width of the angular 
distribution as trajectories with emission angle away from 
the most probable direction start mainly from tha points 
with Z— values away from the mean value. Upon the inclusion 
of the nuclear interaction, such trajectories might 
undergo substantial changes resulting in a modified angular 
distribution. !Ihe fraction of trajectories thus etffected may 
be pretty large in the present case of compact SPG. Phis 
may demand alterations in the SPG obtained by the point 
charge model. Once we incorporate the correct nuclear 
deformations at the scission point, the above mentioned 
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effect can cause substantial changes in the predictions of 
the trajectory calculations based on point charge model even 
with higher D- values. 

¥e have performed trajectory calculations with an 
improved model incorporating the nuclear interaction. The 
calculations include the effects of the fragment deformations. 

The restriction on the Y-coordinate -of the alpha- 
particle (Pig. 5.1) has been removed and all the SPO-parajaeters 
have been allowed a free variation. 

The weight factor used in selecting the SPC from the 

phase space covered by the calculations, has been modified 

to provide a proper time scale with its origin coinciding 

with the instant of scission. Por this purpose the data of 

16 5 

Oheifetz et al. on the -^He emission in fission has been 
used. 

5 .3A Specifying the Initial Oonfignration 

The parameters defining the SPG are the same as used 
in the point charge model calculation (Pig. 5.1). Rotational 
and vibrational degrees of freedom are ignored. The fragment 
deformations at the scission point have been approximated by 
their ground state values. The heavy fragment is taken to be 
spherical as it is having a nearly closed shell while the 
light fragment has been taJcen to be of ellipsoidal shape. 

The deformation of the li^t fragment is specified by the 
ratio of axes along and at right angles to the symmetry axis, 
which is given a value equal to 1,3. 
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^alues of the parameters (initial total kinetic 
energy of the fragments ), 1° (initial kinetic energy of the 
alpha-particle), D (the distance between the centres of mass 
of the fragments) and X (the position of the alpha-particle 
along the fission axis) have been varied over a wide range 
in uniform steps. 'Values of and E^. were chosen randomly 
within narrow ranges around their mean values, which were 
varied in uniform steps. Such a smearing of the values of 
initial kinetic energies is done to provide* an uniform and 
unbiased coverage to the initial phase space of the three 
particles. The parameters Y (the distance of the alpha- 
particle from the fission axis) and (the angle of 
emission of the alpha-particle) are allowed values between 
0.0 fm to 3.6 fm and 45° to 135° respectively. These 
values are chosen randomly for each set of values of para- 
meters Bp, B°., D and Z. 

The range of variation of D is calculated for each 
set of values of B^ and B°, using the constraint due to the 
experimental distribution in the total kinetic energy of 
the three particles. 

As described later, for certain configurations of 
the fragments and alpha-particle, the latter may get trapped 
by the fragments due to the remnant nuclear attraction. In 
particular alpha-particle gets absorbed when the initial 
position of the alpha-particle is very close to one of the 
fragments. Range of the parameter Z has been calculated 
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for different D-values as a function of Y for typical and 
Values by finding the closest distance frora the fragments 
to which the alpha-particle can be placed without getting 
absorbed by the fragments. The range of 1, thus determined, 
is found to be dependent on the interfragment distance, being 
more for larger D values, (However, it might be pointed out 
that this dependence may become milder, once we include the 
correct deformations at the scission point ). 

5 . 3B The Eorces Among the Particles and the Equations of 
Motion 

Apart from the monopole term in the Coulomb force 
among the particles, quadrupole interaction has been included, 
which is derived from the quadrupole interaction potential 
due to a nucleus with ellipsoidal deformation. The potential 
is given in terms of the coordinates (r,9) by the following 
equation ; 

5^(003 e) 

Here c and, a are respectively the values of the semi- 
axes along and at right angles to the symmetry axis of the 
ellipsoid. The symmetry axis is talsen to lie along 6 = 0. 

Ze represents the charge of the nucleus. 

The nuclear force between the alpha-particle and the 
fragments is calculated from the potential obtained by the 
energy density formalism. The nuclear interaction among the 
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17 

fragments has been ignored. The formula given by K'go et al. 


Was used for calculating the nuclear potential V-^-'(S) which 

given below for the spherical nuclei : 


V®) = 

1 2 ^ (3) 

(Ap 

where 

Vs) = 

-Vq exp(-0.27 S^)f S > 0 


= 

-Yq + 6.5 'j S < 0 

with 

O 

11 

27 MeY 

and 

s 

R - ro + A2^^) 

where 

R 

Distance between the centroids of the two 



nucl ei 


""o 

0.98 fin 


I! 

CM 

Atomic masses of the nuclei. 


In case 

of the light fragment, the necessary correc 


tions were made in the expression for Vjj(S) to account for 
the deformation. The coefficient of Uj^(S) in the expression 
for is modified to include the effect of position 

dependent curvature of the nuclsai* surface. For a given 
position of the a— particle, point on the surface of the 

light fragment lying closest to the it is determined (see 
for the method of calculation). If (r,©) are 


appendix 
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the coordinates of this point then the coefficient is given 
as 

. 1/3 , 1/3 

with 

p = L^lisiVt-i pi/3 

p tan"^e + 1 

where p = c/a for the light fragment and refers to its 
mass. 

The value of 'S' is calculated hy using this point 

'P' on the surface of the light fragment so as to get the 

shortest distance between the surfaces of the nuclei. The 

equations of motion are similar to those used in the point 

charge model calculations. The procedure for numerical 

integration is same as described in section 5.20 except that 

-2'3 

the time interval for the steps of integration is made 10 
sec. , moment the nuclear force between the nuclei becomes 
significant. This ensures a good computation accuracy. 

The energy conservation during the entire trajectory was 
better than 1 */. . 

5.30 Procedure of Obtaining the Scission Point Configuration 
(3PC) 

The general method for obtaining the SPC is based on the 
arguments presented in Chapter IV and similar to the one used 
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and a neutron with half life of = 8x10' 


• 22 


in the point charge model calculations. 

1 fi 

The experimental results of Cheifetz et al. on the 

5 

emission of the particle unstable He ane used in the selec- 
tion of the proper time origin ( which coincides with the 
acual scission point). The particle He decays into He 

sec. The 

decay half life is small compared to the time needed for 
the full acceleration of the particles. Thus the decay 
occurs while the light particle still experiences the Coulomb 
field of the fragments. ITom the neutron energy, the enei-gy 

5 

of the particle He at the time of its decay was deduced by 
Cheifetz et al. This provides an independent means to put 
a constraint on the SPO as it demands that the light particle 
should have an energy of 6,3 +0.8 MeV after an average 
time of scission point. As the energy of 

5 

the He at the instant of its decay is essentiiJ-ly calcula- 
ted from the neutron energy, in the present calculations we 

5 

have calculated energy of the neutron emitted from the He 
decay. 


Por each set of the SPC-parameters along with the 
alpha-particle trajectory, a He trajectory was calculated. 
The ^He trajectory was followed for about twice the half life 
of ^He-decay. The energy of the neutron from the ^He was 
calculated by assuming the emission angle to lie within + 5° 
(in the centre of mass frame) at ten equal time intervals. 

An average energy of the neutron was evaluated for the 
exponential decay of ^He. Under the assumption that the SPO 
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4 5 

is identical for the He and He accompanied fission, a 
weight factor W was assigned to the corresponding alpha- 
particle trajectory, which is given as 


¥ = 



.P 


B 


n 


where P^ , P^ , P^, represent the prohahility factors 
corresponding to the experimental distributions in the 
energy of the fragments (Bp), the alp ha- par tide (B^) and 
the neutron from the ^He breakup respectively. 


5. 3D Results and Discussions 

Pigs. 5.12 and 5.13 provide a typical illustration 
of the effects of incorporating the nuclear forces. The 
drastic effects on the final values of the dynaiuical 
variables are found due to the presence of the nuclear 
interaction.. However, this occurs for certain initial 
configurations only. Por some initial configurations, the 
alpha-particle trajectory does not approach the fragments 
close enough and the resulting asymptotic values of the 
dynamical variables are very close to the values obtained 
from the calculations without nuclear forces. In addition, 
some initial configurations result, due to nuclear attrac- 
tion, in the reabsorption of the fragments (see Pig. 5.13). 

In the absence of the nuclear interaction, the very same 
configurations result in the emission of the alpha-particle 
with reasonable asymptotic values for the dynamical variables. 




5.12. The asymptotic angle 0 «l of the alpha particle for two values of 0j_ as a 
* function of coordinate X. Values of D=23,5frn, ES = 1.8MeV, Ep=19MeV» 

Y = 0.0fm. Values of the Ecc, Epore nearly constant and arc given in the 
insets except for the two points on the right hand top corner. •-Corresponds 
to calculations 'without nuclear interaction while +• corresponds to calculations 
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This indicates that the approximation, such as i^pnoring the 
nuclear interaction, erroneously permits inclusion of non- 
physical parts of the initial phase space of the three 
particles in the trajectory calculations. Further illustra- 
tions of the effects on the final values of dynamical varia- 
bles due to inclusion of the nuclear force are provided later 
in the discussions of the tarjectories that result in the 
polar emission of the alpha-particle. 

In the following, we shall be first presenting the 
results obtained by alloT>o.ng all the SPG-parameters to have 


their normal distribution as demanded by the experimental 
Bp and distributions used for assigning the probabi- 
lity weight. The modifications in the results due to inclu- 
sion of ^Lie data in providing the probability weight factor 
have boon discussed. Later the results obtained by imposing 
restrictions on the variation of the interfragment 
distance and tome results on polar emission >rill le 

pres ented . 


5 . '5D- 1 ilesults with Free Variation in all the SPC-Parametgrs 


and Discussion 


( i ) Angular Distribution and its Gorrelations with the Energy 

Ilg. 5.14 gives the calculated angular distribution 
of the alpha-particle by using Bp) weight as well as 

(B , B^, B ) weight. Inclusion of the ‘E^-term‘ in the 
weight factor results in the narrowing down of the angular 
distribution. Both the weights provide satisfactory agreement 
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as regards the width and the mean valtie for the angalar 

distribution of the alpha-particle with the corresponding 

experimental values. However, Ej,, E^) weight provides 

a width in better agreement with the results of the latest 
2 n 

measurement which give a E^HM »17 . 

Eig. 5.15 gives the calculated average energy of 
the alpha-particle as a function of its angle of 

emission 'itie variation of E^ with corresponding 

to both weights show a minimum around the most probable 
angle of emission as is also observed experimentally. However, 
the minimum value oi Eq, obtained by using (E^, E^,, E^) weight 
is 16 MeV, in agreement with the accepted value. While the 
corresponding value as obtained by using (E^, Bj,) weight is 
nearly 1 MeV above the accepted value. This brings out the 
effectiveness of the 'Bj^-term' in the weight factor for 
selecting proper volume in the phase space. Eigure 5.16 
gives the mean value of as a function of B^, It is found 
to be very mildly dependent upon B^ with a small increasing 

tendency >Jith B^. This result agrees well with the experi- 

10 5 

mental observation . However, the Guet’s result^ of 

decreasing angular width with an increase in fragment energy 

(Bj.) could not be reproduced by using (E^, Bp) or (E^, Bp, B^) 

weights (Eig. 5.17). This probably reflects the relatively 

lesser effectiveness of the 'E^-term' in selecting the proper 

values of the dynamical coordinates of the ’fragments as 

compared to those of the alpha-particle as the neutron 

energy depends more sensitively on the energy of light 
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paxticle as compared to the energy of fragments. 

(ii) Dependence of i on E-p and E™ on E 

ct j j g 

Eig. 5.18 displays the calculated dependence of 
average alpha-particle energy (E^) on the fragment kinetic 
energy (Ey). Also shown is the calculated relation between 
Ep and B^. It is seen that the calculated results are not in 
agreement with the experimental ones. It mi^t be recalled 
that the point charge model calculations aHov'ing a free 
variation in all the SPG-parameters also failed to reproduce 
the epcperimental anti- correlation between Ep and E^, 

( iii) Distributions in and E-p 

Pig. 5.19 gives the calculated distribution in the 
energy of the alpha-particle. The width and the mean values 
are in excellent agreement with the experimental values, 
which are used in calculating the weight factor. 

The calculated distribution in the total kinetic 
energy of the fragments is shown in Pig. 5.20. The calculated 
width are found to agree reasonably with the experimental ' 
results while the calculated mean value is about 1 MeV less 
than the experimental value. Some asymmetry is discernible 
in the calculated Bp-distribution. The source for such an 
asymmetry possibly lies in the approximations made in the 
calculations. Por instance, it might be due to larger values 
of the inter fragment distance being erroneously more favoured 
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in comparison to the lower ones. This happens in the present 
calculations as we have approximated the fragment deforma- 
tions corresponding to the scission point tj their ground 
state values. This results in the availability of rela- 
tively higher volume of the initial phase space for the traje- 
ctories (with larger D-values) that result in lower values 
of Bj,. 

The above argument emphasizes the caution to be 
observed in chosing the phase space volume for different 
inter fragment distances. 

(iv) Distributions in the SPO-Parameters 

Big. 5.21 gives the distribution in the initial 
kinetic energy of the alpha-particle. The distribution is 
peaked at about 1.5 MeV with a high energy tail. The 
distribution in the initial fragment kinetic energy is shown 
in Big. 5 , 22 . The distribution is peaked around 28 MeV with 
a tail in the high energy region. The interfragment 
distance (D) is distributed as shown in Big, 5.23. The 
distribution in D is seen to be fairly wide (BWHi 5 fm) and 
centered around 24.6 fm. 

5 . 3D-II R esults with Restricted Variation of Interfragnent 
Distance and Discussion 

We performed calculations with the interfragment 
distance (D) restricted within a range of + 0.5 fm aroxind 
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values varying from 20 fm to 29 fm. The results indicate 
that a fairly reasonable agreement with the experimental 
data can be obtained if the value of I) is restricted around 
24.9 fin. 

( i ) Angular Distribution and Energy Angle Correlations 

Fig, 5.24 gives the calculated angular distribution 
of the a-particle for different ranges of the interfragment 
distance. The calciilated angular distribution for around 

20.4 fm is seen to be narrower (FWHM % 12°) as compared to 
the experimental angular distribution (FWHM % 18°). On the 
other hand, calculations with higher value around 

28.5 fm yield a bit too wide angular distribution (FV/HM % 30°), 
The calculated angular distribution is seen to be of right 
width (FWHlvi 20°) when 'D’ is restricted ai'ound 24.9 foi. 
Unlike the variation in the calculated angular width with 

’D’ , the mean value of the angular distribution is nearly 
unchanged for different 'D* values and agrees with the 
experimental results. 

Fig. 5.25 gives the calculated angular distribution 
for ’D’ restricted around 24.9 fm for two windows of the 
total fragment kinetic energy (Bp). The mean value of the 
angular distribution is nearly constant while the distribu- 
tion becomes narrower for the higher values of the Bp in 
agreement with the experimental results of Guet et iL. . 
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Fig. 5.26 shows the dependence of the average alpha- 
paf'ticle energy as a function of its angle of emission for 
’D' restricted around 24.9 fta. The minima observed in the 
Value of occurs around the most probable angle of 
emission in agreement with the experimental results. In 
addition the minimum value of (= 15 i^ieV) is also in 
agreement with the experimental value, 

( ii ) Energy Distributions and Correlations 

Pig. 5.27 shows the calculated distributions in the 
total fragment energy (Ej,) and the alpha-particle energy 
(E ). The results shown are for ’D’ restricted around 
24.9 fin and are in nice agreement with the experimental 
data used in calculating the weight factor. Pig. 5.28 gives 
the calculated relationship between average alpha-particle 
kinetic energy E^ and the fragment kinetic energy Ej, for 'D* 
restricted around 24.9 fm. The observed fall in the value 
0 f with Bp and its non-linear nature agree reasonably 
with the experimentally observed trend. The dependence of 
average total fragment kinetic energy (Ep) upon '1^* is 
shown in Pig. 5.29. The observed anti-correlation between 

E-r, and E comes out to be -0,5 in reasonable agreement with 

if oc 

12 

the experimental value (-0.3 +0.1) of Mehta et al. 

( iii ) Distribution in the SPO-Parameters 

Pig, 5,30 gives the distributions in the initial 
total fragment kinetic energy (Ep) and alpha-particle 
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energy (B^) for 'D' restricted around 24.9 fia. Distribu- 
tion in E° is pealced around 1.5 MeY with a high energy tail. 
Distribution in shows a peak around 24 MeY with a high 
energy tail, an observation similar to that, seen in Eig. 
5.22, In the present calculations E^ was initially varied 
upto a maximum of 30 MeY which was extended by about 10 MeY 
to check the falling trend in Ep distribution. Further 
extension of this range^may seem desirable by looking at 
the calculated distribution. It might be mentioned here 
that the number of trajectories added in the calculation 
become fairly large when we include higher Ep. This happens 
as higher Bp imply higher ’D’ values, requiring more number 
of points of emission as the distance of the point of 
emission is varied in uniform steps of 0.1 fm. Even in 
calculations for extending the range of Ep from 30 MeY to 
40 MeY, the number of trajectories increased from nearly 
40,000 to 70,000, It is seen that this does not affect the 
main conclusions drawn from the present calculations. 

Pig, 5,31 gives the distribution in the distance of 
the point of emission along the fission axis from the 
centre of the light fragment. It is peaked around 14 fm 
with a full width at half maximum of about 5 fm. Distribu- 
tion in the distance of point of emission from the fission 
axis is found to be nearly uniform. This reflects the 
relative insensitivity of the final variables on the 
coordinate Y. 
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5.3D-III Hesults on the Polar Alpha-Particles and Discussion 

Due to presence of attractive nulcear forces, some of 
the alpha-particle trajectories are bent in the polar 
direction (which is conventionally taken as within 25 ° from 
the fission axis). In pure Ooulomb force calculations, the 
Coulomb repulsion of the fragments prohibits the alpha- 
particle to emerge in the polar direction (refer to Chapter 
II, section 2.2). 

In the present calculations, we have used a fairly 
fine grid to cover the phase space characterising the SPC. 

A significant fraction of the total number of trajectories 
is found to be bent in the polar region. We associate an 
additional term in the weight factor for the equatorial 
alpha-particles to take account of the polar alpha- 
particles. The equatorial term consisted of a product of 
normalised G-aussian functions corresponding to the distri- 
butions of the alpha-particle energy and total fragment 
energy in the equatorial emission. To this equatorial term 
was added the polar term obtained by taking product of 
normalised Gaussian functions corresponding to the distribu- 
tions of the alpha-particle energy and total fragment energy 
in the polar emission. With this weight factor we calculated 
the angular distribution of the alpha-particle over the range 
covering the equatorial and polar directions. Pigure 5.32 
gives the calculated relative yield per unit solid angle as a 
function of the angle of emission alpha— particle 
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allowing all the SPG parameters free variation. The calcu- 
lated probability over a small range of the 9 . was divided 
by the solid angle corresponding to this range to get the 
yield per unit solid angle. The calculated angular distri- 
bution shows all the important features of the experimental 
result on the angular distribution^®*^^. Sharp peak around 
85° corresponds to the equatorial alpha-particles and has a 
PWHM nearly 20°. The rising portions near the angles 0° 
and 180° reproduce the experimental observation for the - 
polar alpha-particles. The polar to equatorial yield ratio 
for the alpha-particles is found to be 2x10 , This number 

is to be compared with the experimental value^" of 5x10 . 

Agreement is fairly satisfactory considering the uncertain- 
ties associated with the calculated values and the experimental 
data as well as the approximate model used for defining the 
SPG. The yield of the polar alpha-particles emitted along 
the light fragment is more than the yield of the ones emitted 
along the heavy fragment in agreement with the experimental 
trend. However, the ratio of these two yields comes out to 
be 1.5 which is small in comparison to the experimental 
value^*^ of 3.8. Calculated mean energies of the polar 
alpha-particle comes out to be about 27.5 MeV and 28,5 HeV 
for emissions along the light and the heavy fragments respec- 
tively. These values are to be compared with the corresponding 
experimental values which are 25,5 i'^eV and 24.8 MeV respectively. 
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The above results indicate that a reasonable agree- 
ment is achieved between the calculated and experimental 
results for the polar emission within the limitations of the 
model used. Important point to be noted is the fact that a 
single set of calculation has succeeded in reproducing 
several important features of both the equatorial, and the 
polar emissions, 

21 

In an earlier attempt, Dahowski et al. concluded 
that the results on polar emission can be reproduced by 
trajectory calculations including nuclear force only if 
certain dubious assumptions are made regarding the nuclear 
force. The present calculations, however, show that the 
nuclear orbitting model of polar emission deserves further 
investigations. The agreement between calculated and experi- 
mental result is remarkable considering the approximations 
in the calculations. The calculations for the polar alpha- 
particles could not be done for different ranges of the 
inter fragment distance due to poor statistics for them. ¥e 
have not studied the effect of incorporating the friction 

term in the nucle«ar interaction. It may be mentioned in 

21 

passing that the results of Dakowski et al. on polar 

omission indicate an improvement in the agreement between 

calculated and experimental results upon inclusion of the 

% 

friction term. 

The nuclear potential cmaployed in the present calcula- 
tions is an approximate representation of the nuclear 
interaction between the alp ha-par tide and the fragments. 
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However, the main difference between the potential used in 

the present calculations and the optical model potential 

21 

used by Dahowski et al, appears to be associated with 

presehce of the repulsive core in the former. The presence 

of hard-core permits a possibility of a-p article being 

re-emitted during the evolution of the trajectory even if 

the separation between the alpha-particle and a fragment 

goes below the sum of their half density radii. ¥e consider 

alpha-particle to be absorbed only if once re-emitted during 

the evolution of the trajectory it comes back again under 

the attractive part of the nuclear potential to see the 

repulsive core in it (thus the fragment- alp ha configuration 

will resemble a bound state). On the other hand, for some 

polar trajectories, the alpha-particle does not see at all 

the repulsive core of the nuclear potential. The latter 

type of trajectories will resemble the trajectories as 

21 

obtained in the calculations of Dakowski et al, . Pigs, 

5-33 and 5.34 show the trajectories corresponding to the 
above mentioned possibilities. The effect of repulsive core 
in the nuclear potential is manifested in the region of kink 
(marked as HO in Pig. 5.33) in alpha-particle trajectory. 

The experimental results on the alpha-particle scattering at 
energies below Coulomb barrier are not reliably usuable to 
probe the deeper parts of the nuclear potential. The calcula- 
ted results on scattering at such energies are rather insen- 
sitive to these parts of the nuclear potential. Thus the 
experimental data do not provide a clear picture regarding 
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the presence or absence of the repulsive part in the nuclear 
potential. Theoretical investigations for understanding the 
behaviour of the deeper part of nuclear potential suggest 
that repulsive core may arise from the antisymmetrization 
effect etc.^^^^^. 

( i ) On the Validity of Delayed Emission Hypothesis 

Possibility of explaining the polar emission by 

20 

delayed tripartion is discussed by Piasecki et al. . They 
discard this possibility as their calculations indicate 
necessity of a delay of the order of 10“ sec, (implying , 
an increase in inter fragment distance of 200 fm) for 
explaining polar emission. However, considering the alpha- 
particle trajectory shown in iig, 5,55, we can see that if 
the alpha-particle were emitted in the appropriate direction 
around the region marked as ’RG’ , it can come out in the 
polar direction. This will imply a delayed emission at a 
point when fragments are nearly 27 fm apart starting with 
inter fragment distance of ^ 20 fm. Thus we see that the 
alpha-particle trajectories which are bent in the polar 
direction after experiencing the repulsive core of the nuclear 
potential will remain valid trajectories with reasonable 
Values for E and E-^ even after the elimination of the 
repulsive core provided we assume the alpha-particle 
emission to be after a little delay (which may cause an 
increase of about 10 fm in the interfragment distance). 
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It is seen that the presence of nuclear attraction, which 
is ignored in the calculations of Piasecki et al,^^, is 
responsible for bending of such trajectories in the polar 
direction. 

5 • 3E Surnmary and Conclusions 

Ihe present calculations have indicated that substan- 
tial changes in the SPG may be incumbent once we include 
the nuclear interaction among the three particles in the 

light charged particle . accompanied fission. This result 

25 

is in agreement with the findings by Dakowski et al. while 

26 

results of Guet et al. are not consistent with it. 

The calculations with a free variation of all the 
SPC-parameters and including nuclear force have succeeded 
in reproducing quite satisfactorily the angular distribution 
of the alpha-particles, its energy distribution, its average 
energy as a, function of the angle of emission etc. These 
calculations predict some characteristics of the polar 
emission in reasonable agreement with the experimental 
findings. In addition, the calculations predict reasonably 
well the results for the ^He data a* (calculated mean 
neutron energy comes out to be in excellent agreement with 
the experimental value). The results indicate that the 
method used by us which permits a free variation in all the 
SPC-parameters is capable of providing accurate and reliable 
information regarding the SPG, Heed of using a realistic 
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model for describing the SPG is brought out by the present 
calculation. Inclusion of the nuclear forces is one step 
towards making the point charge model more realistic and the 
results obtained are encouraging. Use of ^He-data has come 
out to be fruitful in proper selection of the SPG from the 
phase space of the three particles. Furthermore it improves 
the agreement between the calculated and experimental ' 
results. However, it is seen that the selection of the 
proper phase space by the weight factor obtained by using 

5 

He-data is more effectively done for the alpha-particle 
dynamical coordinates in comparison to that for the fragment 
dynamical coordinates. This is expected as the energy of 

5 

the neutron from He breakup is comparatively less sensitive 
to the fragment motion. Similar type of data, which have 
sensitive dependenbe upon fragment motion (such as, on the 
scission neutronsi), will be quite useful in fixing a proper 
time origin for the trajectory calculations in a more 
reliable way. 

The calculations with no restriction on 'D' failed 
to reproduce observed anticorrelation between iij, and 
and Guet's result of decreasing angular width with increase 
in Ep. These results along with the energy and angular 
distributions and energy-angle correlations could be repro- 
duced by calculations with restricted variation in interfragment 
distance (D). The calculations indicate that best agreement 
with experimental data is achieved when U is restricted 
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axotmd 24.9 fta in a narrow. range. Thus, the condition of 
reproducing the anti correlation and the G-uet’s result is 
seen to necessitate a compact scission configuration (with 
small spread in the value of 1), albeit with large initial 
kinetic energies compared to the result of point charge 
model calculations. 

The present calculations have not included the 
friction term in the nuclear potential. This term may have 
important effects in the results, e.g, , those on the polar 
emission. Also, the initial values of fra^ent deformations 
are taken approximately while their time dependence has not 
been considered. These approximations have affected the 
results of the present calculations. In particular, this 
may mean that the value of 'D’ obtained by calculations 
with restricted D-variation should be taken only tentatively. 
As stated earlier the present calculations are exploratory 
in nature and we find that there is a need of more realistic 
model for describing the SPG. For example, the approxima- 
tions regarding the fragment deformations have resulted in 
an unrealistic bias for higher 'D' values. This is very 
clearly reflected in the assnnmstry observed in the calculated 
distribution of the total fragment kinetic energy with no 
restriction on D (Fig. 5.20). . The mean value is also 

found to be lower by a small amount ( 1 MeV) from the 
experimental value. This bias for higher values appears 
due to comparatively higher range available for the variation 
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of the X-ooordinate (see Pig, 5.1) of the alpha-pahticle. 

The bias can be reduced if we constrain the calculations to 
the trajectories which have value of in a narrow range. 

A calculation with restricted between 80° and 85° shows 
that the mean value of the Bj, becomes in excellent agreement 
with the experimental value and also the asymmetry is substan- 
tially reduced. A calculation with these approximations 
removed will be free from such unrealistic bias coming in 
the selection of the phase space of the three particles to 
characterise the SPO. 
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CHAPTER - VI 


■ SIR'MARY AED COHCLUSIOiJS 

The phencraerjon of mission of light charged particles 
(LCP ) is a relatively less understood aspect of nuclear fission 
and deserves further investigation. The present thesis reports 

(i) a ccmparative experimental study of two aspects of the 
li^t charged particle accompanied fission, namely, the 
equatorial and the polar emissions, 

(ii) a study of the trajectory calculations to extract 
information on the dynamical aspects of the fission 
process uti.lizing the characteristics 'of the equatorial 
LCPs and to probe the possibility oi' an unified description 
of the polar and equatorial emissions, of LCPs. 

Short critical surveys of the LCP accompanied fission I 

and the trajectory calculations along with some suggestions and | 
comments are also Included. , 

Simultaneous laeasui’eraents on the equatorial and polar 
LCPs emitted in the iission of hy thermal and fast ( 600 , 

K.eV ) neutrons have been performed using a novel ionization 
chamber geometry, 'ui entification of LCP being accomplished with 
the help of a AH-ii semiconductor detector telescope. The | 

geometry permits ui-e of broad fission source without losing the ; 
angular resolution necessary for an appropriate identification , 
of the polar and equatorial LCPs. Detailed Monte Carlo calcula- | 
tions on the detection efficiency have been performed. Results ; 
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on the absolute and relative yields of the equatorial and polar 
LCPs in the thermal neutron induced fission are in agreement 
with the earlier work. The yields of the polar and eq^uatorial 
protons show a marked increase when the neutron energy changes 
from the thermal to the fast value while the yields of the polar 
and equatorial alpha-particles indicate a mild decrease. Most 
interesting finding of the present study is the observed 
similarity in the change of yields of ICl ;jith a change in neutron 
energy for the equatorial and the polar emissions. Such a result 
provides a rather clear evidence suggesting a similarity in the 
emission mechanisms of the equatorial and polar LCPs in fission. 
Two hypotheses on the mechanism of polar emission, namely, the 
nuclear orbitting hypothesis and the quantum mechanical diffrac- 
tion hypothesis, appear most promising in view of the present 
experiment as both of them attempt essentially to provide an 
unified description o± the polar and equatorial emissions of 
LCPs. 

It might be worthwhile to mention the rare nature of the 
phenomenon that has been studied in the present experiment as we 
are dealing with coincidence measurgnents on neutron induced 
reactions with extromely low cross-sections (for fast neutrons, 
the cross-section is as low as few tens of microbarns for the 

pclar LCPs). 

Traijectorj' cal.culations have been performed under point 
charge model approximation exploiting seme recent results on the 
energy-angle correlation in the alpha-particle accompanied fissior 
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The calculations were* done by providing a restricted vai’iation 
in the interfragrnent distance around values lying between the 
limits expected on the basis of different fission theories. 

Other parameters cliorac terising the' scission point configuration 
were allowed their normal statistical distribution as denanded 
by the experimental weight factors -used for selecting the phase 
volme characterising; the scission point corn igurat ion ■ from the 
initial phase space of the alpha-particle and fragments. A 
relatively compact scission point configuration with interfrag- 
ment distance belo-'r 21 fm could reproduce various energy-angle 
correlations of the alpha-particle and iragnents arid is thus 
favoured by the point charge model calculations. 

An improvement in the trajectory calculations has been 
attempted by 

(i) making tho point charge model more realistic by including 
the effects .,luo to fragement def orraauti ons and the nuclear 
interaction between the alpha-par tide and fragments 

5 ' 

(ii) using the data on the breakup of unstable He particle 
produced in LCP accompanied fission /.or fixing a proper 
time origin in the calculations co-inciding with the 
scission moMieiit 

(iii) using comparatively finer gride for various initial 
variables characterising the scission point configuration 
providing' oji uniform coverage of the Initial phase space 
of the three particles and looking ior a possible unified 
description c:C the eq_uatorial and pold omissions of ICPs. 
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With these improvements, calculations have been performed 
with a free variation of all the initial variables characterising 
the scission point configuration; proper selection from the phase 
space being done by the probability weight factor obtained by 
using the exp erimentcil energy distributions of the alpha-particle, 

n: 

the Iragements and the neutrons from the '"^He breakup. The 
results obtained reproduced very well the angular distribution 
of the alpha-particle, the distributions in the energy of the 

alpha-particl e and fragaents and seme energy-angle correlations. 

5 

Use of data cna the He - breakup is found to be useful in getting 
over the basic uncertainty regarding the proper fixing of the 
time origin in trajectory calculations but need of more accurate 
data is felt. Hovrever, the observed anticorrelation between 
energies of the alpha-particle and fragments and the observed 
narrowing of the alpha-particle angular distribution with 
increasing fragment energy could be reproduced only when a 
restriction in the variation of the interfragment distance is 
' put. Calculations with restricted variations in interfragment 
distance have been performed covering low as well as hi^ values 
of interfragment diiStf-nce. These calculations indicate that a 
reasonable agreement between the experimental and calculated 
cnergy-angLe correlations and distribution is possible when the 
Interfragment distance is restricted around 24.9 fm. 

Thus it is rGa].ised that significant changes in the 
reconstructed scisnton point configuration may be incumbent once; 
the point charge model is made more realistic by incorporating 
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nuclear interaction etc. Need of further improvement has been 
suggested by a consideration of some of the calculated results. 
They concern with inclusion of accurate deformations at the 
scission point as well as their time dependence during the time 
evolution of the system. In addition friction term in the 
nuclear potential etc. should be considered. It is felt that 
incorporation of these improvements might affect the conclusions 
drawn from the present calculations. 

Reasonable reproduction of the experimental characteristics 
of the polar emission of alpha-particle is achieved by the present 
trajectory calculations with improved model. 




fig. A.1. Schematic diagram showir>g the alpha -parti 
and the light fragment used in calculating the 
shortest distance between their equivalent 
sharp surfaces . 
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Bqn. (4) is a fourth degree alg'ebraic equation and has been 
solved numerically on the computer, riewton-.iaphson method 
has been used for this purpose. The initijii guone for the 
value of k is obtained by approximating point 'i ’ with the 
point D suown in the Pig. This value k is tr,iven as 

kg = [/Tx^ + ^ 

where ’r^' is the nuclear radius parameter and *A’ is the 
mass number of the light fragment. 




